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This  report  constitutes  the  fifth  part  of  an  article  on 
Water  Waves  being  prepared  for  the  new  edition  of  the  Encyclo 
paedia  of  Physics  (Handbuch  der  Physik)  published  by  Springer 
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F.  EXACT  SOLUTIONS 


F 


The  word  "exact  in  this  context  is  generally  understood  to 
mean  solutions  in  which  there  has  been  no  approximation  in  the 
equations  or  boundary  conditions.  However,  this  usage  of  the 
word  does  not  exclude  neglect  of  viscosity  and,  in  fact,  since 
positive  results  have  been  obtained  only  for  perfect  fluids,  the 
treatment  below  will  be  restricted  to  them.  Indeed,  the  present 
results  in  the  theory  of  exact  solutions  are  restricted,  with 
few  exceptions,  to  a  very  special  class  of  motions,  namely,  those 
which  can  be  represented  as  steady  two-dimensional  flows. 

In  section  32  some  general  theorems  will  be  established.  In 
section  33  waves  of  maximum  ampl i tude-to- length  ratio  are  dis¬ 
cussed;  because  the  methods  are  intimately  related,  we  have  also 
included  in  this  section  a  discussion  of  Havelock's  method  of 
approximating  periodic  waves,  Section  34  treats  methods  of  ob¬ 
taining  explicit  exact  solutions  and  of  various  ones  which  have 
been  obtained.  In  section  35,  the  last,  existence  theorems  are 
discussed,  but  only  in  a  descriptive  way,  for  proofs  are  highly 
technical  and  lengthy. 
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#  32,  Soma  general,  theorems 

This  section  will  be  devoted  t.o  several  theorems  of  a  rather 
general  nature  concerning  the  motion  of  a  fluid  with  free  surface 
in  a  gravitational  field.  The  theorems  In  32 «  are  mostly  of  a 
kinematics,  1  nature  and  are  associated  with  the  phenomenon  of 
mass  transport  already  discussed  in  section  27 qc  .  The  last 
part  of  this  section  is  devoted  to  several,  theorems  on  energy 
and  momentum.  In  section  32 /j  some  theorems  concerning  waves  in 
heterogeneous  fluids  will  be  established.  In  33y . .  several  dif¬ 

ferent  ways  of  formulating  the  problem  of  motion  with  a  free  sur¬ 
face  will  be  described. 

3  2  at.  Kinematical  theorems --mass  transport — energy  integrals 

The  first  theorem,  due  to  M.  S.  Longue t-Higg ins  [1953],  Is 

n 

independent  of  the  presence  of  a  free  surface  or  of  the  nature  of 
the  time  dependence.  Let  {())  =  'f + <  f'  describe  a  space-periodic 
motion,  i.e.  -*  ffj.)  .  The  definition  of  tf>  will  be  nor¬ 

malized  so  that 


O,  },tjcx  =  o 


<o  (32.  1) 
Note  that  if  this  condition  holds  for  on©  value  of  y  ,  If  holds 
for  all  since 
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In  equation  (2.10")  we  shall  write 
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(32.2) 
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In  the  following  we  define  an  average  by 


(32.3) 


Theorem:  In  a  non-uniform  space-periodic  motion  u^r1, ;  -  j» 

each  decrease  with  increasing  depth  ,  provided  either  c p(x-Ut)=Q 

or  i'""  4*  =  O  . 


This  may  be  proved  as  follows.  Consider  first 


uv-t  vv 


Then 
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By  a  similar  computation  it  follows  that 


2>  «>• 
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f1  =  f  [  (V  ■*  >o' 


(32.4) 


(32.5) 


since  we  have  assumed  that  is  not  constant.  It  is  evident 

from  (32.4)  that,  if  the  fluid  is  bounded  below  by  ^  s  -  i,  ,  then 


(32  6) 


if  it  is  infinitely  deep,  it  is  an  assumed  bounds  - y  condition 
that  cf>  ->o  as  £-*-coand  hence 
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h  V  as  j -’-“o 


(32.7) 


In  either  case  it  then  follows  from  (32.5)  that  ft ^  is  an 

Increasing  function  of  ^  and  hence 


3  _ 

K  f  '  °> 


(32.8) 


with  equality  occurring  only  for  fact,  even  more  caa 

be  concluded,  for  (32.5)  is  like  itself  with  cp  replaced  by 
2cpx  .  Hence,  by  repeating  the  above  reasoning  one  may  establish 


that 
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(32.9) 


Next  consider  uv  -  V1  .  A  similar  computation  s!-  jws  that 


**  '  >  ]0 


Hence 


-  V1  =  C(t)  -  ^  |  , 

or 


(32. 10) 


It  follows  from  (32.8)  that 


[  f  4  C]  ^  Vvs  ~[  f  -C)  ,  '  (*j 
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(32.11) 


are  each  increasing  functions  of  ^  ,  i.e.,  they  decrease  with 


increasing  depth.  For  infinite  depth  Longuet-Higgins  shows 


further  that,  if  axes  are  chosen  such  that  w*  o  at  ^s~*f 


then 


F 


th«  qun.ntl.ties  I  m'|  ,  |  t/i  and  ,|?4|  all  decrease  exponentially  to 
*ero.  He  had  shown  earlier  ( 1950]  for  exact  waves  (we  shall  not 
carry  through  the  proof)  that 

i.  _  ',7^ . 

2 

(3  2,  12) 

Hence  it  follows  that, 

-T  |  -  -  JL.  '7^ 

U  ;  }  5-k  "  I  3tv  *  *  (32.  13) 

For  purely  progressiva  waves  this  quantity  vanishes,  but  we  re¬ 
call  that  for  standing  waves  we  found  earlier  a  constant  pressure 
fluctuation  of  double  the  wave  frequency  [see  (27.62)  and  (27.65)] 
if  second-order  terms  were  retained. 

lass  transport.  In  section  27  [see  (27.39)  and  (27.41) ]  it 
was  shown  that  a  forward  drift,  called  "mass  transport”,  occurred 
in  progressive  waves  if  second-order  terms  were  taken  into 
account.  It  waa  shown  by  Rayleigh  [1.876]  in  a  proof  valid  only 
for  infinitely  deep  iltiid  that  mass  transport  must  always  occur. 
The  proof  :1s  independent  of  the  dynamical  free -surface  condition. 
Levl-Civita  [1912]  and  later  Ursa  11  [1953]  developed!  methods  of 
analysis  to  include  both  finite  depth  and  nonperiodic  waves; 
essentially  Ursell's  analysis  has  also  been  given  by  Nekrasov 
[1951]  for  infinite  depth.  The  analysis  given  below  is  due  to 
Longue t -Higgins  [1953]  and  is  similar  to  that  used  in  the  pre¬ 
ceding  theorem.  We  note  that  Starr  [1945]  has  also  given  an 
instructive  and  perspicuous  derivation  of  Rayleigh's  theorem 
for  infinite  depth.. 
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Take  the  wave  as  moving  to  the  left  with  velocity  c  (in  the 


sense  of  section  7)  and  inpose  a  uniform  velocity  c  in  the  oppo 


site  direction,  so  that  the  motion  is  reduced  to  a  steady  one 


generally  in  the  positive  x -direction  in  the  sense  that  w > £  »  o 


We  m a;,  then  write  the  complex  potential  in  the  form 
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where  F*€  >  t  >  0 
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We  take  ^ 
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at  a  crest  and  assume  y  *0 
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(32.  16) 


Denote  by  T(^)  the  time  required  for  a  given  particle  to  progress 


f 


In  the  original 


wave  motion,  the  time  elapsed  between  the  passage  of  two  sue 


cessive  crests  over  a  given  point  is  X/c 


ielt< 


>  X  /  c 


the 


particle  is  being  transported  with  the  wave  and  it  will  be 


reasonable  to  call 


u(t) 
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(32. 17) 


th®  mass  transport  In  the  direction  of  wav®  notion.  The  follow 


ing  theorem  is  true 
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Theorem :  Both  T  and  U  decrease  with  increasing  depth,  and, 
with  the  assumed  definition  of  £  ,  U  >0  . 

The  theorem  may  be  proved  by  the  following  computation: 


fi(\)  < 

Tf  t)  =  I,  f/' 
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2s 
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(32.  18) 
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***** 


(32.  19) 
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(32. 20) 


The  details  of  the  computation  are  almost  identical  with  those 
used  in  deriving  (32.4)  and  (32.5).  Since 


\  "i* ^  1  0 

it  follows  from  (32.19)  that  Then,  since  T'(i^J>0 

unless  the  flow  is  uniform,  it  follows  that 


In 
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I' 


(32.  21) 


T‘(r )  *  O, 

with  equality  holding  only  for  .  As  in  the  earlier  theorem, 

the  computations  can  be  exterded  to  yield 


W)>0,  ? 


It  now  follows  immediately  from  (32.17)  that 

u'(  *')  £  0  t 


lliiiilll!iii|lllli"i  III !  > 


(32.  22) 
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(32.23) 

. . . 


with  the  equality  holding  only  for  the  bottom  at re aw line.  If 
the  fluid  la  Infinitely  dee#,  then  U'  >0  for  alt  .  To  com¬ 
plete  the  proof  we  must  show  that  U  >0,  If  the  fluid  Is  in¬ 
finitely  deep,  it  its  evident  that 


lm  T(f]=i'  02.20 

Hence  )jm  U  =0  as  and  the  conclusion  follows  from 

If  the  depth  is  finite,  we  compute 


"T  f*  Q) 


,L, 


(32.25) 


here  use  has  been  made  of  the  Schwarx-Bunyakovakii  inequality. 

(We  have  written  >  rather  than  ^ ,  for  the  equal  sign  will  hold 
only  in  the  trivial  case  of  a  uniform  flew. )  1  now  follows  that 

UM)|>0  and  hence  that 


U  >  O 


(32.  26) 


since  0  O  .  This  completes  the  proof  of  the  theorem,. 

The  method  of  analysis  can  foe  extended  to  prove  an  analogous 
theorem  for  nonperiodic  steady  motions  which  approach  uniform 
flows  as  x  -*  ±  °o  ,  in  particular,  to  the  solitary  wave. 


Momentum  and  energy  integrals.  We  close  this  section  with 
several  momentum  and  energy  integrals,  most  of  which  have  been 
found  by  Levi-Civita  [1912,  1921),  Starr  [ 1947ab,  1948)  and  Starr 
and  Platsaan  [1948]. 

Let,  us  again  take  the  wave  as  moving  to  the  left  without 
change  of  form  and  impose  an  opposite  velocity  o  which  brings 


the  profile  to  rest  (or,  equivalently,  consider  the  motion  re¬ 
lative  to  a  coordinate  System  moving  with  the  wave).  I*#1©  t  >  t  h  © 
velocity  potential  be  as  in  (32.14).  Consider  the  area  bounded 
by  two  streamlines  and  ~  ,  say  £  ■*.  y,  (x)  and 

and  two  vertical  lines  a  wave  length  X  apart.  To  this  area 
apply  the  theorem 
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(32.27) 


This  yields 
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since  Cf>n  -  O  on  the  streamlines.  Moreover,  since  Cp( X  t ), }  l) 
c)  ♦  'v'x,  })  ,  =  c  ■+ w  and  y.  (x4 ))  *  7; (x)  ,  the 


right-hand  side  of  (32.28)  may  be  written  as 
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Expanding  and  rearranging  give 
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If  one  now  applies  the  operator  A  <lx  to  (32,30) 


(32.30) 


one 


obtains 
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or,  after  multiplying  by  y  f  and  rearranging, 
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(32,32) 


i.e.,the  kinetic  energy  per  wave  length  between  two  streamlines 
equals  j  c  times  the  momentum  in  the  direction  of  the  wave 
(here  to  the  left). 

Next  let  us  write  the  integral  (2. 10')  in  the  form 


t  f  yf  +  |*  -  f 


(32.33) 


the  fcrm  of  the  constant  having  been  chosen  for  later  convenience 
Write  the  terms  t  -r  f  *  >  as  follows  : 

-  -V-V+  ~  (32.34) 

Here  we  have  used  the  second  equation  of  (2.6).  We  may  now  write 
(32.33)  as  follows 


(32.35) 


Next  let  us  integrate  equation  (32.35)  over  the  same  area  as  is 
described  In  the  preceding  paragraph.  First  consider  Pf'jv-J/Dt 
Since  the  motion  is  steady  in  the  selected  coordinate  system, 
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where  the  last  equality  follows  from  the  continuity  equation. 
Hence 


f[  =  ^v-fu4t^).Q  dj  - 


A  l  - 


(32.36) 


since  ^  sO  on  the  streamline  boundaries  and  the  Integrals  over 
the  vertical  boundaries  cancel  from  periodicity.  The  integrated 


equation  then  becomes 


pfuW*)* i?  -t  C  J]  ?u  jcr 
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(32.37) 


I  f  one  eliminates  the  second  integral  by  means  of  (32.32),  one 


obtains 


rX  :  pxv 


r  A 

iitr  -t  i  K  j  J* 
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(32.38) 


Equation  (32.38)  has  a  simpler  aspect  if  the  two  streamlines 
are  taken  as  the  free  surface  7  (V)  and  the  bottom  --j  t -■  lv  .  Then 
f) ( x ,  =  0  and  the  third  integral  becomes 

,  >  ,•> 
j  P}?VM  i*  +  i  j,  l*  . 


Moreover. 


c  * 


l*  -  0 


(32.39) 


if  the  x-a*is  is  taken  at  the  mean  water  level.  This  follows 


from  the  following  sequence  of  equations,  similar  to  thou*  used 


in  (32.36): 


t  tr 


(32.  40) 


Equation  (32.39)  now  allows  us  to  give  a,  simple  physical  inter¬ 
pretation  of  the  constant  cf  in  (32.33) .  For  if  (32.33)  is  in¬ 
tegrated  along  '-i  i  -  h  ,  and  account  is  taken  of  (32.39),  one 


finds 


I 

)\ 


Cw* 


(32.41) 


■<L 

i.e.  C,  is  the  mean  square  velocity  of  fluid  along  the  bottom. 
The  inequality  follows  easily  from 


|  z  '  '  tfJjr-l  )lx  =  “O* 

6  ’  (32,42) 

If  the  fluid  is  infinitely  deep,  w  -*  O  as  -*  -  eo  ,  and  (32.41) 
reduces  to 

V.  *W 

C  “  Cf  '  (32.43) 

If,  following  (15.27),  we  let  ,  ]fAV  ,  . 

*  ^  *’ 

Tfld„  denote  the  average  kinetic  energy,  the  contributions  to 

this  due  to  the  x  and  y  velocity  components,  potential  energy, 

and  momentum  in  the  direction,  of  wave  motion,  respectively,  then 

# 
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(32.32)  and  (32.88)  may  be  expressed  as  follows: 


2  y  -  t'rn 

A* 


=214. 


(32.44) 


where  the  last  term  of  the  second  equation  is  0  for  li  =  no  . 

The  first  equation  is  essentially  due  to  Levi-Civita  [1912,  1921] 
the  second  to  Starr  [ 1947b], 

W©  note  another  simple  consequence  of  (32. 41) ,  due  to  Levi- 
Civita  [1924].  Let  us  integrate  (32.33)  along  the  free  surface 
for  a  wavelength  and  divide  by  ~  p \  .  Then,  remembering  our 

choice  of  x-axis  as  the  mean  water  level,  we  find 


m  “  •  i 


(32.45) 


On  the  other  hand,  if  we  compute  the  velocity  at  the  intersection 
of  the  mean  water  level  and  the  profile,  we  also  find 


(32. 46) 


Hence  the  absolute  value  of  the  velocity  at  the  mean  water  level 
equals  the  root-mean-square  velocity  along  the  surface  profile 
or  along  the  bottom,  or,  indeed,  along  any  streamline,  for  in 
the  reasoning  in  (32.40)  we  could  have  substituted  any  stream¬ 
line  ^ for  "jjs-li  and  obtained 


(32.47) 


Starr  and  Platzman  [1948]  have  used  the  relations  above  to 
derive  some  general  relations  concerning  the  flow  of  energy  in  a 
periodic  wave,  le  recall  that  the  average  flux  of  energy  in  the 


direction  of  wav*  motion  Is  given  by  [cf.  §8,  (IS. S3)  and  (15.87)] 
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ft  follows  from  the  second  formula  In  (32,44)  that 


hi* 


3  Tj,,  ~  i  x  '*■ tV)  ^ 

ience,  with  £a*  :  ,  we  obtain  from  (32.48) 


<32. 49) 


_L  . 

4 


i  £  r,. 

I  l  f. 

w  £ 


-a* 


(32. 50) 


'This  should  be  compared  with  the  result  derived  in  section 


Iot  infinitesimal  waves  with  neglect  of  surface  tension  [cf 

/ 


(15.25)  and  (15.26)],  namely,  :tt  =  g  £*„  .  Equation  (32.50)  is 
consistent  with  this,  for  to  the  order  of  approximation  involved, 
iJy-  4*  »nd  C,  st  .  However,  for  waves  of  finite  height  it 
was  shown  in  section  27  at  [cf,  eqs,  (27.42),  (27.43)]  that  to 

■i*  • 

the  second  order  of  approximation  4*  >  K*  .  Plataman  f 1947  ]  has 
verified  that  this  remains  true  when  4 th-order  terms  are  kept. 

Several  of  the  above  results  have  analogues  for  steady 
motion  of  nonperiodic  waves,  provided  that  n(*)*+c?  as  x  — *t* 

In  such  a  way  that  ,  V x  1®  finite.  Under  such  circumstances 


V  m  %» 

C -  a  C  and  the  following  results  may  be  established  (the  nota- 


tion  is  that  of  (15.31)  with  obvious  extensions): 
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(32.51) 
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For  details  of  the  proof  one  may  refer  to  Starr  [1947b].  From 
the  last  two  equations  follows 

£V  -  =tfU  “»*  i  4.T.,  Ut  >  ?*•  (32.52) 

We  note  that  the  second  equation  of  (32.51)  is  a  special  case  of 
a  more  general  one  applying  to  any  steady  motion: 


ix)  -  'x  ( x)  -  *r*  -  tj-it. 

»  ■) 


(32. 53) 


where  the  constant  is  zero  under  the  conditions  of  (32.51).  The 
proof  is  analogous  to  that  of  (8.6).  Here 


7M 

f  -t 


j  etc. 


32/rt,  .  Waves  in  heterogeneous  fluids 

The  first  two  theorems  proved  below  are  also  true  for 

# 

homogeneous  fluids  and  were  first  proved  for  this  case.  The 
last  theorems  deal  specifically  with  heterogeneous  fluids.  In 
the  extended  form  they  are  all  due  to  Dubre i 1-Jacotln  [1932]. 

A  flow  will  be  called  barotropic  if  both  the  pressure  and 

f) 

density  are  constant  along  streamlines.  We  first  derive  the 
energy  integral  for  such  flows.  The  equations  (2.6)  may  be  writ¬ 
ten  in  the  following  form  in  two  dimensions: 
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Since  J>  is  assumed  constant;  on  a  streamline,  <*f>„  +v"Pj 
follows  from  (32.27)  and  the  definition  of  £  that 
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(32. 55) 


In  particular,  if  the  flow  is  steady^  b  is  also  constant  along  a 
streamline.  For  steady  flow  it  is  a  consequence  of  the  incom¬ 
pressibility  condition  that  f  is  also  constant  along  a  stream- 
1  ine . 

The  following  theorem  was  proved  by  Burnside  [1915]  for  a 
homogeneous  fluid.  He  gives  two  proofs,  of  which  the  second  can 

be  carried  over  to  the  present  more  general  situation  with  no 
change.  It  will  perhaps  give  more  substance  to  the  theorem  if 
we  remark  that  Gerstner's  wave  (see  section  34/J  ),  which  is  not 
irrotational,  satisfies  the  other  conditions  of  the  theorem. 

Theorem :  The  only  steady  two-dimensional  irrotational  mo¬ 

tion  of  a  fluid  subject  to  gravity  for  which  all  streamlines  are 
also  lines  of  constant  pressure  is  a  uniform  flow. 

Let  the  streamlines  be  given  by  z  const.  Since,  from 

the  remark  following  (32.55),  E  -  const,  along  a  streamline, 
we  may  write 


(32 . 56) 


(Burnside  shows  that  one  may  generalize  (32.56)  by  replacing  <j  j 
by  a  function  j,(y) .)  Since  the  motion  is  irrotational,  A  V'-O 
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But  then 


A  Jdjj  [  Sty)  "  =  C', 


which  yields  after  some  computation 


2  ^  t  '(f)  “  i  !  £  "  J  $ )  [  £  '  “  '  c“  *  3 } )  f  "  J  4  J 


We  write  this  in  the  form 


fu  C'":  >)  "  C\  (  1 


It  then  follows  from  (32.56)  and  (32.58)  that 


'8.,  "  ^  '  6  .  f-jj.  ~  ^  ”•  o 
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But  then 
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(32. 57) 


(32,  58) 


and  is  a  function  of  ^  only.  Hence,  since  A  ^  -  VC ,  s  0  ,  ^ 


is  a  constant  and  the  flow  is  uniform. 


The  next  theorem  was  first  proved  by  Levi-Civita  [1925]  for 
homogeneous  fluids.  Fenchel  [1931]  showed  that  his  hypotheses 
could  be  weakened  and  Dubrei  i-Jacotin  [1932]  extended  Fenchel 's 


proof  to  heterogeneous  fluids.  The  gist;  of  the  theorem  is  that 


I 


if  the  surface  profile  moves  withe  at  change  of  form,  then  the 
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whole  velocity  field  is  steady  In  a  coordinate  system  moving  with 
the  surface.  The  theorem  will  be  formulated,  in  the  moving  co¬ 
ordinate  system. 


Theorem ;  Let  a  possibly  heterogeneous  fluid,  bounded  below 
by  a  horizontal  plane  ^  -  b  ,  be  flowing  irrotal  iona  1  ly  in  the 

x-direction  with  discharge  rate  Q  it)  and  with  a  fixed  surface 
profile  ^  -  *7»>)  .  If  and  »  satisfy  the  conditions 


U|  <  *5  <  4 


>  £  >  0. 


1  '  J  '  (32.50) 

then  the  velocity  potential  f  Jj)  is  independent  of  t  . 

First  we  derive  a  boundary  condition  at  the  free  surface. 
From  the  condition  of  constant  pressure  and  the  assumption  that 
the  surface  profile  is  an  invariant  streamline  it  follows  that 


U  -+-  lr  =  O  on  ‘'l'  Z  Q 

e>  > 

It  then  follows  as  in  (32.55)  that 


-  qy  4  JPjiL  sr  0  on  1 Ps-o  , 

Dt  "  d  D  t  (32.60) 


IZowfver,  this  conclusion  holds  now  only  on  this  one  streamline. 

The  complex  potential  maps  the  region  of  the 

z -plane  occupied  by  fluid  onto  the  strip  0  ,  where  the 

free  surface  corresponds  to  '/  r  o  ,  the  bottom  to  -  - Q  and 

X  = -±  oo  to  (p  =  -£  e©  .  Let  F(Y);  '£ 4 >  "fy  be  the  mapping,  unique  up 
to  an  additive  real  constant,  of  the  fluid  region  onto  the  strip 
-I  sf/  &  O  with  /:  i«  corresponding  to  sr  foo  „  Then 


iiliilllllillHiiliiHihHii . . 


||  II  l|n|l  III  II  ill  ll  'I1  II . I' 11  1 


<9(hF(y) 


Iff  I  Hi 


(32.61) 


1 . . . . . . . . . I . . 


j 

I 


I  H 

f-m 


evidently  satisfies  the  requirements  for  1 a»d*  in  ,act 
is  determined  uniquely,  up  to  the  added  constant  In  F  »  *3> 
and  rj(»)  .  Now  substitute  into  (32 , 6°) 
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(32, 62) 


which  we  may  write  In  the  form 

Q*  4  A  Q  ■+  3  -S  ^  =  C  (32.63) 

where  A  and  3  are  independent  of  *  .  Division  by  4%  4  ^ 
is  possible  since  (32.59)  Implies  that  this  does  not  vanish.  Note 

also  that 


3 


»  < 


(32. 64) 


and  that  both  A  and  3  may  be  considered  as  functions  of  $  • 
Consider  two  cases:  (a)  A  -  const.,  <b)  A  *  const.  (a)  In  thi 
case,  since  3  is  independent  of  t  and  Q  4  A  Q  is  independent 


Of  <$  ,  it  follows  from  (32.03)  that  both  equal  constants.  It 
now  follows  from  (32.64)  that,  unless  this  constant  is  zero,  the 


profile  y{x\  will  be  unbounded  and  the  first  part  of  (32.59)  will 
be  contradicted.  Hence,  in  case  (a)  *)  -  const,  and  the  mapping 
F  must  be  of  the  form  -*  =  -t  •£  ,  *  and  A  real.  It  then  .toll  on. > 
that  A  and  hence  Q*  -  O  ,  i.e.  the  flow  is  uniform.  (b)  Let 
,  A."  be  two  different  values  of  A  ,  At  •  Write  equa¬ 

tion  (32.63)  for  each  value  and  subtract.  This  yields 


3,  -  3v 
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(32.65) 


But  then  Q  is  evidently  independent  of  t  ,  Hence  also 
H},t)  '=  Q  F(yj  is  also  independent  of  t  .  This  completes  the 
proof. 

The  next  theorem,  due  to  Dubrei t- Jacotin  [1932],  specifi¬ 
cally  requires  that  the  fluid  be  heterogeneous. 

IllfifiCflU :  Suppose  the  motion  of  an  Incompressible  hetero¬ 
geneous  fluid  to  be  irrotat ional ,  the  free  surface  to  move  with¬ 
out  change  of  form,  and  that,  in  a  coordinate  system  moving  with 
the  surface,  conditions  (32.59)  are  satisfied.  Then  not  only  is 
the  velocity  field  steady,  but  also  £  ,  p  and  P  are  constant 
a  1  ong  t  he  s  t.  re  am  1  i  nes  . 

It  follows  from  the  preceding  theorem,  that  the  velocity  is 
steady,  hence  that  £  =  However,  we  may  still  concei . ably 

have  j-j-s.  p \X'  ;  £,  ,  |»  * ?  (x  t  The  equations  (32.54)  may  now 

be  written  in  the  form 
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two  equations 
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ie  assume  that  the  corresponding  functional  relations  may  be 
solved  and  write 


where,  from  (32.66) 


(32. 67) 
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From  the  equation  expressing  lncompressibi lily,  namely 
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(32, 69) 


We  shall  assume  l)f  /d£  $  0  everywhere,  and  may  thus  write 


<3  (E,  Y)  _  pt  f£,f) 

i/^2)  P*  (7jJ 


(32.70) 


Since  the  left-hand  side  is  independent  of  t  ,  it  follows  from 


the  form  of  the  right-hand  side  that  we  may  set  both  sides  equal 

1 0  A  f  £  )  »  i .  e  . 
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(32.71) 


Let  us  suppose  h'&)$  0  ,  e.g.  ^(£fj  2  0  ■  Then  j 9  must  be  a  func¬ 

tion  of  the  form 
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(32. 72) 


in  some  neighborhood  of  £f  .  If  h  c)  vanishes  for  some  values 
of  E  ,  let  f  be  the  first  zero  larger  than  Et  ■  Then  from 
(32.71)  and  (32.72) 


Jct  h{t)i 
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(32.73) 


But  (32,73)  can  hold  for  all  t  only  if  f'  -O  ,  i.e.  if 
p  const. j which  is  contrary  to  the  assumed  heterogeneity. 

Moreover,  at  least  one  such  zero  of  exists,  for  we  already 


know  from  (,<12,60)  that  £  is  constant  along  the  free?  surface, 

so  that  in  steady  motion  the  Jacobian  in  (32.70)  vanishes  for 

if' *  0  .  Hence  4Vc)  =  0  for  the  corresponding  value  of  E  , 

W©  must  conclude  that  A  ' E 2  Q  .  This  implies,  from  (32.70)  that 

£  -  fi  l  i/)  and  p  s  f,=j  •  from  (32.68)  and  the  condition  t  z.  Q 

1 1 

on  the  free  surface,  ft  follows  that  also  p*  .  Hence  , 

f  and  £  are  all  constant  on  streamlines. 

The  last  in  this  complex  of  theorems  is  also  due  to  Dub* . . 

Jacotin  [1932] 

Theorem :  There  cannot  exist  ir rotational  waves  in  a  hetero¬ 
geneous  fluid  such  that  the  profile  is  propagated  without  change 
of  form, 

Thls  . . .  lows  Immediately  from  the  first  and  last  theorems 

proved  above,  and  is,  of  course,  subject  to  condition  (32.59). 

Ihls  striking  result  is  a  1 1  the  more  s-  * n  view  of  the  fact  that 
Gerstner’s  wave  (section  34/J)  does  .  de  a  steadily  propa¬ 
gating  wave,  even  in  a  heterogeneous  f.uid.  The  theorem  also 
casts  some  doubt  upon  the  significance  of  the  linearized  theory 
of  irrotational  wave  motion  in  a  heterogeneous  fluid  as  developed, 
for  example,  in  Lamb  [1932,  |  235],  Such  a  wave  evidently  cannot: 
be  considered  as  a  first  approximation  to  an  exact  steady  solution. 

32  t  •  S^e=trans forma t i ons . of  the  boundary-value  prob lero 

By  means  of  introduction  of  new  variables  or  other  devices, 
it  is  possible  to  formulate  the  boundary-value  problem  for  exact 
solutions  in  a  variety  of  ways.  Several  such  formulations  will 
be  considered  in  section  3  4  exon  inverse  methods.  Here  we  give  a, 
few  which  seem  to  be  of  general  interest. 


Inversion  of  ffjr)  .  One  elementary  but  important  transfor¬ 
mation  has  already  been  introduced  in  section  32«-in  the  discus¬ 
sion  of  mass  transport.  This  is  the  inversion  of  the  velocity 
potential  f (}/)  when  >f'\  vanishes  nowhere  within  the  fluid 
and  treatment  of  +  as  the  independent  variable.  This  has  the 
advantage  that  under  certain  circumstances  the  domain  of  defi¬ 
nition  of  the  independent  variable  can  be  given  exactly:  when 
^  is  the  Independent  variable,  the  domain  of  definition  is 
one  of  the  unknowns  of  the  problem.  For  example,  if  the  motion 
is  reducible  to  a  steady  flow  with  discharge  rate  Q  ,  one  may 
take  the  surface  profile  to  correspond  to  ,'A'=0  and  the  bottom 
streamline  to  correspond  to  -  —  Q  Hence  the  domain  of 

definition  of  \ '  *■  Is  the  strip  0±Y  '  '5  .  if  the  fluid 

V 

is  infinitely  deep,  the  domain  is  the  hai  lane  V'  -  0  .  When¬ 
ever  !  can  be  taken  as  the  independent  variable,  then  one  can 

ft  p 

also  express  *  +•  as  a  function  of  T  .  It  has  been  estab¬ 
lished  independently  by  Gerber  [1951]  and  Lewy  [1952a]  that  the 
equation  describing  the  free  surface,  £  -  ^  )  ;  is  an 

analytic  function  of  '-f>  at  all  points  for  which  i >s4  0  ■ 

Stokes'  "second  method".  In  the  introduction  to  section  27 
it  was  mentiosied  that  Stokes  [  1880],  in  a  supplement  to  an  ear¬ 
lier  paper  in  his  collected  works,  developed  a  method  for  approxi¬ 
mating  exact  periodic  waves  which  is  different  from  the  straight¬ 
forward  generalization  of  in  f ini tesimal -wave  theory  expounded  in 
that  section.  This  method  is  based  upon  use  of  f  as  the  inde¬ 
pendent  variable  and  expansion  of  3  as  a  Fourier  series  in  -f  : 
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for  Infinite  and  finite  depth  respectively;  the  «„  may  be  taken 
to  be  real.  Here  ^  is  taken  as  in  the  preceding  paragraph.  The 
coefficients  A*  are  tc  be  determined  from  the  condition  that 
the  pressure  be  constant  on  the  surface,  i.e.  from 
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If  the  mean  water  level  is  taken  at  u  ■=.  q  and  the  fluid  is  infi¬ 
nitely  deep,  then  C  :  i‘  ;  we  shall  consider  on ly  this  case  here, 
Then  equation  (J2.7G)  may  be  expressed  as 
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(32. 77) 


Substitution  of  (32.74)  in  (32.77)  yields 
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After  multiplying  the  two  factors  and  reducing  the  cosine  products 
to  cosines  of  sums  and  differences,  the  resulting  expression  may 
be  put  into  the  for* 
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where  the  -4„,  i  and  Cv are  forms  of  the  third  degree  in  the  am  \  , 
The  coefficients  of  the  individual  cosine  terms  must  then  be 

equated  to  zero.  This  results  in  an  infinite  sequence  of  equa . 

t ions ,  each  involving  all  the  a.,  S  and  -j.  ■'  ;  “ -  •  In  order  to 

proceed  further,  one  must  devise  some  method  for  approximate  de¬ 
termination  of  the  A,  :,s  .  Stokes’  procedure  was  to  assume  that 
each  a„  could  be  expanded  In  a  power  series  in  some  parameter, 
the  initial  term  in  the  series  having  the  power  t  .  This  allows 
one  to  carry  through  a  step-by-step  improvement  in  the  approxima¬ 
tion  of  the  Xr  ‘S  by  including  successively  higher  powers  oi  the 
parameter.  We  shall  not  pursue  the  matter  further,  but  remark 
that  the  most  systematic  arrangement  of  such  computations  seems 
to  have  been  devised  by  Sretenskii  [1952]. 


Le vi -Ci vita’s  dif f e rential- di ffe rence  equation.  T h e  f o 1 1 o w - 
lug  theorem,  due  to  Levi-Civita  [1907],  reduces  determination  ot 
u/iff)  for  steady  flow  over  a  horizontal  bottom  to  solution  of  a 
differential  -di.  f  fe  rence  equation  . 

Theorem :  The  complex  velocity  z.  j-»  v  of  an  irrotationa  l 

gravity  flow  with  constant  discharge  rate  3  and  with  u  »  f  > O 
must  satisfy  the  differential-difference  equation 
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Conversely,  anj  function  -v  f)  satisfying  (32.80)  which  is  regu¬ 
lar  in  the  strip  -i'Sslft  m2,  finite  at  x  ,  real  on  Jm  f 3  '  1 
and  has  a  >  t  ■>  o  represents  such  a  flow. 

In  order  to  derive  (32. 80) ,  we  note  first  that  the  functions 
im  !  + )  and  '■  '  ’  *  *  ►  both  have  vanishing  real  parts  for  y~  -  1*  and 
consequent!)  can  be  extended  b>  reflection  to  the  strip  •*■!  /  i\* 
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The  f ree-sur face  condiL.on  maj  be  expressed  b>  the  equation 
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or,  by  making  use  of  the  extended  definitions  oi  »'  and 
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Consider  the  '  -nc  1 1 or 
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Evidently,  H  Is  defined  and  is  regular  on  the  line  V 


and 


thus  in  some  neighborhood  of  this  line.  From  (32.83)  it  follows 

if*  i  _ 


that  H  (  J  -  )  rC)  ,  hence  that.  H  f)  =  0  in.  its  .region  of  defi¬ 

nition.  Equation  (32.80)  follows  from  the  fact  that  j/ !  ■£  r  <  C4)  ~ 

1/ w (f  1  <  O')  .  For  proof  of  the  converse  we  refer  to  Levi  -Clvita ’s 

paper.  Levi-Clvlta  also  gives  a  special  form  of  (32.80)  appro¬ 
priate  to  a  space-periodic  flow.  Cisotti  [1919]  generalized  the 
preceding  theorem  to  include  a  variable  discharge  rate.  The 
equation  (32.80)  may  be  considered  to  contain  the  equation  (22.30), 
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when  in  that  equation  ,  in  the  sense  that  linear!- 

i at  ion  of  (32.80)  by  assuming 


„a/'  C  '  t 


yields  (22.30). 


Rudzki * s  transformation.  The  following  transformation  was 


a ppm re n 


tly  first  introduced  by  Rudzki  (1898).  It  has  later  been 


used  by  many  others  in  the  investigation  of  exact  water  waves. 

The  validity  of  the  reformulated  boundary  condition  is  not  limited 
to  periodic  waves.  However,  it  is  assumed  that  a  coordinate 
system  has  been  selected  with  respect  to  which  the  flow  is  steady. 
It  is  again  assumed  that  u  >  £  >  o  .  Let  9  be  the  angle  between 
the  velocity  vector  aid  the  positive  '.-axis.  Then  one 

may  write 
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where 
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(32. 86) 


Here  c  is  some  typical  velocity,  say  the  wave  velocity  as  defined 
in  section  7.  We  consider  o  as  a  function  of  {  and  let  /'  =0 
correspond  to  the  free  surface.  The  free-surface  condition  may 
then  be  expressed  by 
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and,  from  (32  86) 
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Hence  (32.87)  becomes 
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or,  since  from  the  Cauchy-Rlemann  equations, 
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It  one  can  find  a  function  io  v  regular  in  the  strip  a&lfr 
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with  I'SI'-J  £  >  and  with  its  real  and  imaginary  parts  satisfy 


ln«  (32  88)  or  on  ‘^*0  .  one  may  then  construct  from 
it  a  free-surface  flow  with  gravity.  Of  course,  further  con¬ 
ditions  must  be  imposed  at  Vs  -  Q  or  as  ’/  -*  ~ao  . 


—  ras  •  The  following  transformation  is 

due  to  Nekrasov  [1921,  1951].  It  will  be  assumed  that  the  sur¬ 
face  is  periodic  with  period  X  ,  symmetric  about  a  crest  and 


that  the  fluid  is  infinitely  deep  and  Jim  ur 

/  -*  -  «£ 


Let  the 

origin  in  the  z-plane  Lc  taken  at  a  crest,  Q  be  the  free 
surface ,  and  assume  u>f  >0  .  In  addition  to  the  z-  and  f- 
p lanes,  we  introduce  a  T -plane, 


C--  —  pe 

related  to  the  f-plane  through 
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With  a  cut  along  the  negative  £  -axis*  there  is  a  one-to-one  cor¬ 
respondence'  between  the  various  don, a  ins  CAOBD  shown  in  Figure  46. 


FIGURE  46 


The  relation  between  the  z-  and  5 -planes  will  be  deter¬ 
mined  by 
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where  '*  l  f]  is  regular  in  the  disc  and  is  related  to 


by 
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(32.  93) 


The  form  of  b  shown  in  (32.92)  follows  from  the  assumed  proper¬ 
ties  of  the  motion.  Since  j5  =  l  on  the  free  surface,  the  con¬ 
dition  of  constant  pressure  may  be  expressed  by 
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and  satisfying  (32.96). 


From  such  a  function  on©  can  easily  con¬ 


struct;  a 
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flow  with  free  surface. 


rasov's  Integral,  equation.  Nekrasov  also  considers  the 
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function  of  (32.92),  but  as  a  function  of 
from  (32.88)  and  compute 
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One  may  formally  Integrate  this  equation  and  obtain 
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where  I  is  the  Integration  constant;  yK  is  related  to  the 


velocity  at  the  crest, 
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Substitution  of  <32.98)  into  (32.99)  yields  the  following  ©qua 


tion  for  the  relation  between  v  and  #  on  tie  boundary 
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(It  follows  from  (32.98)  that  the  denominator  does  not  vanish). 
It  is  known  from  the  theory  of  functions  of  a  complex  variable 
[see,  e.g.,  Cara  the odory ,  Funkt  ionentheorie,  Bd .  1,  §  1.47-9, 

Birkhauser,  Basel,  1930]  that,  if  a  function  is  regular  within 
and  on  a  closed  Jordan  curve,  it  is  determined  up  to  an  addi¬ 
tive  constant  by  giving  either  its  real  or  imaginary  part  on  the 
boundary.  In  particular,  in  the  case  at  hand  we  may  express  the 
value  of  9  on  the  boundary  (  z  /  in  terms  of  r  on  the  boundary 
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12. 102) 

From  the  assumed  symmetry  about  a  crest  follows  t%4)  -  -  ?  Y/l )  »  so 
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that  (32.10°'  may  be  expressed  as  follows; 
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Substitution  of (32. 100)  into  (32.103)  yields  Nekrasov's  nonlinear 
integral  equation  for  &(/)  : 
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(32.  104) 
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e. 


K  one  can  llnd  Q  satisfying  (32. 104) ,  one  can  thou  reconstruct 
l.  and  hence  the  whole  flow. 


Nekrasov  [1.928,  1951 J  carried  through  a  similar  analysis 

when  the  depth  Is  finite.  We  shall  only  sketch  It.  In  Figure 

35  suppose  that  £=  represents  the  bottom  (  *  is  not  the 

mean  depth)  and  if*' =  -  eg  the  corresponding  streamline.  In  the 
T 

3  -Plane  this  maps  into  a  circle  of  radius  fQ<i  ,  where 


Po  -  a  ■ 

In  (32.92)  Uf)  becomes  a  Laurent  series.  The 
for  Q(f)  remains  the  same  in  form  as  (32.104), 
function  log,*  !  is  now  replaced  by 


(32 . 105) 


integral  equation 
but  the  kernel 
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Moiseev  [ 1957b]  has  further  generalized  Nekrasov's  equation  so 
as  to  allow  a.  wavy  bottom. 

The  solution  9(  {)  of  (32.104)  will,  of  course,  depend 
upon  the  parameter  ^  ,  except  for  the  trivial,  solution  9  SO 
corresponding  to  a  uniform  flow.  It  is  possible  to  show  that 
not  all  are  allowable.  Let 
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It  then  follows  from  (32.102)  that 
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hence  that 
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From  this  follow'* 
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V1 1 1. 1  a  t  ’a  lot  e  g  r  a  1  e  q  u  ail  o  n .  Even  though  we  shall  not  consi¬ 
der  Its  contents  In  any  detail,  it  would  be  improper  not  to  men¬ 
tion  an  important  paper  of  VI llat  [1915].  Vi llat  wished  to  find 
the  steady  motion  of  a  fluid  in  a  canal  of  given  bottom  profile 


and  also  with  a  given  top  profile  over  the  part  of  the  fluid  up¬ 


stream  of  some  point.  Downstream  of  this  point  the  top  profile 
is  one  of  constant  pressure .  The  boundary  condition  on  the 
free  surface,  (32.89),  is  modified  by  introduction  of  new  vari- 
ableSj  and  a  pair  of  integral  equations,  one  of  them  nonlinear, 
is  derived.  The  method  is  also  applicable  if  the  upstream  "cover" 
is  absent;  and,  in  fact,  becomes  a  little  simpler.  The  chief  use 
made  of  the  procedure  by  Villat  is  as  an  inverse  method  in  which 
the  free  surface  Is  given  and  the  corresponding  bottom  and  cover 
determined. 

3 3 .  Waves  of^  maximum  amplitude 

In  the  higher-order  theory  of  infinitesimal  waves  one  of  the 
important  effects  of  including  higher-order  terms  was  to  make  the 
profile  more  peaked  at  the  crests  and  flatter  in  the  troughs. 

The  effect  was  the  same  for  either  steady  progressive  waves  or 
standing  waves.  Since  the  peakedness  increased  with  increase  of 
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the  miipl i tude-to-wave- length  ratio,  It  seems  reasonable  to  con¬ 
jecture  that  there  is  some  bound  to  this  ratio  and  that,  if  a 
wave  01  ampl  i  tude- to- length  ratio  exists,  it  will  be 

characterized  by  a  corner  or  a  cusp  at  the  crest,  at  least;  if 
capillarity  is  neglected.  it  has  never  been  proved  that  such 
waves  exist.  However,  if  one  assumes  their  existence,  it;  is 

possible  to  prove  some  necessary  properties.  This  will,  be  done 
be  low. 

Following  an  earlier  erroneous  investigation  of  Rank  me 
(  IS 65 J ,  Stokes  I IfiflO,  p.  225]  showed  that,  if  a  corner  occurs 
in  steady  motion,  the  angle  included  between  the  tangents  must 
be  1*0  .  Michel  1  [  1893]  assumed  that  a  periodic  highest  pro- 

. . .  we  cxist*  *nd  showed  how  to  compute  the  coefficients 

of  an  associated  series,  but  without  proving  convergence, 

Havelock  1,1919]  made  Michel  1  ’s  procedure  the  basis  of  a  general 
method  of  approximation  to  periodic  progressive  waves,  again 
with  no  proof  of  convergence.  Micheil’s  wave  was  later  investi¬ 
gated  by  a  different  procedure  by  Nekrasov  [1920).  However, 

Nekrasov  did  not  carry  his  computations  to  the  same  degree  of 
accuracy  as  Michel  1  and  Havelock,  so  that  the  numerical  results 
are  discrepant.  More  recently  Yamada  [1957]  rediscovered 
Nekrasov’s  method  and  carried  through  the  calculations  with  the 
necessary  accuracy;  the  results  are  now  in  substantial  agreement 
with  those  of  Havelock  and  Michel  1. 

Penney  and  Price  [1952b],  in  their  work  on  standing  waves 
of  finite  amplitude,  include  an  analysis  intended  to  show  that, 
if  there  exists  a  standing  wave  of  maximum  amplitude  with  a 
corner  at  the  crest,  then  the  angle  must  be  90°.  G. I,  Taylor[l953 ] 
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has  questioned  the  validity  of  the  proof,  and  it  appears,  In  fact, 
to  be  untenable.  On  the  other  hand,  in  the  sane  paper  Taylor 
reports  the  results  of  experiments  which  appear  to  confirm  Penny 
and  Price's  prediction.  in  view  of  the  present  unsatisfactory 
state  of  the  theory,  it  will  not  be  further  discussed  here. 


Stokes'  theorem.  We  prove  firs*-  Stokes'  theorem  on  the 
angle  at  a  corner  in  steady  flow.  Let  the  corner  be  at  the  orl- 
R*n  J  “  0  >  free  surface  be  the  streamline  -  ©  ,  and  -/■'  -  O 

at  the  corner.  Since  "]  zzc  is  assumed  to  be  a  corner,  it  must 
also  be  a  stagnation  point  and  the  constant-pressure  condition 
on  the  surface  may  be  taken  in  the  form 


?-  '  ° '  (33.1) 

In  the  mapping  from  the  y  to  the  f-piane  the  point  —  o  must 

be  a  branch  point,  so  that  in  the  neighborhood  of  •>  r  o  the  com- 

& 

plex  velocity  potential  will  take  the  form 

/  .  v. 

~  \  '  (33.2) 

If  <  O  is  the  angle  between  the  right-hand  tangent  to  the 
corner  and  OX,  then  near  ?  =0  equation  (33.1)  can  be  written 
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This  can  hold  for  all  small  r  only  if 
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It  also  follows  that,  if  Q(_  is  the  angle  between  the  left-hand 
tangent  and  OX,  then  Sieoy.  and  oc.  =  fifO*  ~  <a<  +  so  that  the 

surface  Is  symmetrical  about  OX  near  the  corner.  If  ) k  %  0  cor¬ 
responds  to  the  region  occupied  by  fluid  and  if  the  branch  of 
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with  - J* if  *, at <:  ^  -  is  taken,  then  the  complex  velocity 

0  : 


ht ' 

potential  has  the  following  form  near  •: 

JT 


.*<**>  *  -  a  i  ) 

3  *  r  *  *1 


2/x 


JA  * 


> 

i*-y) 


*  S,  *• 


( 33 . 4 ) 


The  streamline  =  0  has  a  corner  at  J  =-0  with  included  angle 

Q 

120  ,  In  this  case  the  flow  is  to  the  right.  The  inversion  of 


(33.4)  gives 


for  f  near  0. 


VJ 


-<7rA  vj 

e  -f 


(33.5) 


33  c/ 


Period  Lc 


•  aye  of  maximum  height 


Let  us  suppose  that  a  periodic  progressive  wave  of  maximum 


amplitude-length  ratio  exists.  We  may  take  this  as  a  steady  flow 
with  complex  velocity  potential  +  ^  ~  Cp  ti  F  and  with 


(33.6) 

Let  the  origin  of  the  z-plane  be  at  one  of  the  crests,  the  sur¬ 
face  profile  correspond  to  ']{  -  O  ,  and  the  origin  of  the  f-plane 
t0  that  oi  the  z-plane.  Then  the  free  surface  condition  may  be 
taken  in  the  form  (33.1). 


I'rr,  \'f-A  -  c. 
j  — b  -  so 


method.  First  we  give  Michel  l 's  procedure  for 
finding  ■  As  we  have  done  earlier,  we  shall  write 


1* 


c 

c 


(33.7) 
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From  the  assumed  periodicity  and  symmetry,  fy  is  an  odd  periodic 
function  of  (f>  with  period  c,\  for  /'O  ,  From  (33.7)  follows 

T  *\  . 

*> '  ’  “f  J  -  -  .—7  ♦  2  2i- 

•i**  ■  <  (33.8) 

For  f  O  ,  is  an  even  periodic  function  of  #  with  re¬ 

movable  singularities  at  the  crests:  we  expand  it  in  a  Fourier 
series : 


2)  0 

2)<p 


Tf  \  _ 

f  \  I  ^  0  +  •'?! 

w*  /I  L* 


cl 

c  A 
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<33. 9) 


The  <Xk  are  real.  Substitute  (33.9)  into  (33.8)  and  rewrite  it 
in  the  following  way: 


r  =  j 


-j 


Sir 


2  r  ~  Cp 
c  A 


Now  consider  the  function 


(33. 10) 


-  I  is  -n 


/  C  A 


(33.  11) 

Z(f)  is  defined  in  the  whole  lower  half-plane,  is  regular  for 
^  <  O  ,  and,  as  J(--> -co  ,  2  I  f )  ~>-»naa/c,\.  Moreover,  from  (33.10) 
Z  is  also  real  on  the  real  axis  and  hence  may  be  extended  by  re¬ 
flection  to  the  upper  half-plane.  Z  Is  then  a  function  with  only- 


singularities  on  the  real  axis  at  the  points  <£***>  associated 
with  the  crests.  The  form  of  the  singularity  may  be  determined 
from  (33.5) .  In  fact,  near  f-O 


II 


7T  W 
dr  5 


—  Jirr  --  p 

'if 


(33.12) 


Hence  2(f)  has  singularities  of  the  form 


”  J 

f  -  r  t 

,  /I  -  o,  1 1,  x  2.,  •  *  f 

(33. 13) 

along 

the  real  axis, 

and  only  these,  so 

that  it  must 

have  the 

form 

Z( f)  --  - 

— ---  cot 

j  c  ,\ 
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(33.  M) 

From 

(33.  14)  2 

t: / ic  /'  as  Z  -  '-'jo 

Then  from 

( 33 .11) 
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Since 

Z  must:  be  real 

for  real 

f-  ,  it  follows  that 

- 
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3 

(33. 15) 

and 
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TT 

3c> 

- 

cA 

(33.16) 

If  now  follows  from  the  definition  of  Z(f)  that 


77  + 
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•  i  Air  f/c 
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(33.17) 

which  yields,  after  integration,  inversion  of  the  logarithm  and 
use  of  (33.6)  to  evaluate  a  multiplicative  constant, 


v'Ct)  = 

<■’  c  Vi 


|  jf 


y  Sin 
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r . 3*? 


wb#rt  /,  s  I  and  the  ir  are  real.  The  branch  of  the  root  . . . 


be  chosen  so  that  its  argument  lies  between  for  ^ ~  C 

From  (33.  li)  one  finds  immediately  also. 


rfprc^Te 
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33.  19) 


Aside  frown  the  first  one,  the  coefficients  In  (33.18)  or 


are  still  to  be  determined.  The  constant-pressure  condi¬ 


tion  for  the  surface  profile  is  still  available  for  this  purpose, 


for  we  have  made  use  of  the  equation  (33.4)  or  (33.5)  only  through 


the  value  of  the  exponent.  The  value  of  the  gravitation  constant 


has  not  entered  into  (33.18)  or  (33.19).  In  fact,  a  comparison 


of  (33.5)  after  differentiation  and  (33.18)  in  the  neighborhood 


o  f  -f 


yields  immediately 


£  I  -t  6 ,  t  ^  *  1  j  — 


1  4  C  tCv 


(33.20) 


so  that,  once  the  4j-r,  or  Cr  are  determined,  the  relation  between 
wave  length  and  velocity  may  be  found.  This  method  could  pre¬ 
sumably  be  pursued  to  obtain  a  sequence  of  further  equations  for 
determination  of  the  .  However,  Michell  proceeds  sonewMat 

differently.  If  we  differentiate  (33.1)  with  respect  to  $  ,  we 
may  write  the  free  surface  condition  as  follows  (cf.  (32.54)  and 


following)  : 


JL  a 
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it  11  for  f  =  a 

F  5r 


U  f'  for  t-0 


(33.21) 
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Substitution  .of  (33.19)  in  (33.21)  yields  an  equation  of 


the  following  form 
4 


1 

J 


7 7  l 


,w*  r4  [ 

A  c  A  / 


—  -f  /4,10s  l-Z$  4 
c.  A  1  e  A 


I* 


(33.22) 


where  the  depend  linearly  upon  the  C*'s  ,  'and  the  <4„  depend 
upon  them  in  a  more  complicated  manner.  The  derivation  of  (33,22), 
especially  of  the  right-hand  part,  and  of  the  particular  dependence 
of  the  A„  *  and  3*  S  upon  the  c.* 's  is  rather  tedious  and  we  refer 
to  either  lie hell's  original  paper  or  preferably  to  Havelock's 
more  general  and  systematic  treatment.  Equating  coefficients  of 
the  individual  cosine  terms  leads  to  a  set  of  equations  relating 
£.v/<j  )t  ,  C,  ,  Cx  ,  •  •  • .  The  values  as  computed  by  Havelock, 
which  we  assume  to  be  somewhat  more  accurate  than  Michel l * s  own, 
are  as  follows: 

Q  ) 

L~  z  0.1 11- In  ,  C  -0.04M;  ct-  0.0(14  tc3  0  o®mtC4~  0,0014. 

cl  '  (33.23) 

The  value  for  ^  )i  icx  should  be  compared  with  that  for  infinitesi¬ 
mal  waves,  namely  2  rr  .  Substitution  of  j  c  ^  for  f  in  (33.  in) 
yields  the  velocity  at  a  trough: 


U  —  C  VT  [  )  *  C(  *  C-v-Cj-*  J  —  2l9  c- .  (33.24) 

From  +  2^  so  one  may  now  find  for  the  trough  and  hence 
the  amplitude -wave  length  ratio: 


(33.25) 


f-  H  I 


H.  Jeffreys  [1951]  has  recently  reexamined  the  basis  of  the 
Michell-Havolock  method  of  approximation  and  concludes  that  an 
apparert  discrepancy  between  the  values  in  (33,23)  and  equation 
(33.20)  does  not  really  indicate  a  numerical  error  in  the  compu¬ 
tations.  * 

We  note  in  passing  that  Michel  1  also  gave  the  form  of  i  '£) 
analogous  to  (33.19)  which  must  hold  if  a  highest  wave  with  cor¬ 
ner  exists  in  a  fluid  of  finite  depth. 


Method  of  Nekrasov  and  Yamada.  This  method  makes  use  of  the 
^  -plane  introduced  In  (32.57)  and  related  to  the  f -plane  by 
(32.58).  We  may  again  make  use  of  Figure  46  but  must  keep  in 
Blind  that  In  the  z-plane  there  is  now  a  corner  at  0  with  an  in¬ 
cluded  angle  of  1 2.0° .  Hence  (32.59),  the  equation  relating  the 


and  y -planes ,  must  be  replaced  by 


Jl  wo 


(33. 26) 


and  (32.60)  by 


u r 


if  --  , 

^3" 


Q-i)  3 
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(33.27) 

The  coefficients  are  now  to  be  determined  by  the  con¬ 
stant-pressure  condition  on  the  free  surface  taken  in  the  form 


(33.28) 
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one  obtains  as  the  equation  analogous  to  (32,63) 
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(33,29) 


(33.30) 


This  yields  a  set  of  equations  for  determination  of  ^  /i  /cv  , 

A(  » Aj,  i.  *  *  *  .  The  actual  computation  appears  to  be  as  tedious 
as  that  of  Michel l's  method  and,  in  fact,  Nekrasov's  [ 1920)  com¬ 
putations  do  not  seem  to  have  yielded  as  accurate  results  as 
Ml che 11 's.  However,  as  mentioned  earlier,  Yamada  [1957]  has  set 
up  a  systematic  computation  procedure  and  has  obtained  results 
in  substantial  agreement  with  those  of  Michell  and  Havelock. 

Once  and  the  have  been  determined,  the  surface  profile 

can  be  found  in  parametric  form  by  integrating  (33.29)  with  res¬ 
pect  to  /  from  0  to  '/  .  Figure  47,  reproduced  from  Yamada ' s 
cited  paper,  shows  the  form  of  the  profile. 
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FIGURE  47 
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33(1  .  Have  Lock’s  approximation  for  gravity  mvob. 

In  a  paper  already  cited  several  tinea  above  Havelock  [1,919] 
extended  Michel L's  method  of  construction  of  periodic  waves  of 
naxinun  amplitude,  outlined  in  the  preceding  section,  to  one  for 
construction  of  periodic  waves  of  any  allowable  amplitude -length 
ratio.  Up  to  the  present,  no  one  has  proved  the  series  involved 
to  converge.  However,  as  Havelock  points  out,  the  method  has 
attractive  theoretical  features:  the  parameter  describing  the 
family  of  waves  occurs  in  the  form  e  "  where  /l  varies  from  0, 
corresponding  to  the  highest  wave,  to  oo  ,  corresponding  to  in¬ 
finitesimal  waves. 

The  method  starts  out  exactly  1  ike  Michel l's  up  to  equation 

(33.19)  except  that  it  Is  not  assumed  that  =  O  corresponds  to 
the  free  surface.  We  recall  that  in  Michel,  l’s  analysis  the  con¬ 
stant-pressure  condition  did  not  enter  completely  until  after 

(33.19) ,  in  particular,  in  (33.21).  Havelock  assumes  instead 

that  this  condition  is  to  be  satisfied  on  some  other  streamline, 
'if'  ss  -  at  ,  which  will  then  be  taken  to  correspond  to  the  free  sur¬ 
face.  The  condition  may  still  be  written  In  the  form  (33.21) 
provided  that  one  replaces  'Y  -  O  by  'j/  ~  ~  ^  ,  For  ^  -  cc  one 


may  write 


«:5 


(33.31) 


,  -ItJTctf/fA 

where  Q  ,  the  being  the  same  as  those  in 

(33.19).  Havelock  shows  that  one  way  express  c  -  /} i  in  the 

following  form 


. ""Hi . . . .  "*"■ 


(33,32) 


and 


It*  I  in  the  form 


1,1 
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s , ,,  z_?  r  s .  +  s,„  'll  ‘  fg  £  os  Zt  +  £ 

t''  '  c).  c  A  J  1  '  ti 


3-# 

■Oj  f  o  s  — -»  ■ 
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(33.33) 


Here  the  A„  5  are  rather  coo.pl icated  expressions  in  the  H e  .5 
but  also  involve  cosh  ~~>  /  c  ■  linearly;  the  3,  1  Rre  linear  ex . 

press ions  in  the  i  A  with  coefficients  which  are  functions  of 

-l"Ol/c)t  ,  ,  ...  Q  * 

e  .  Havelock  finds  complete  expressions  for  the  S  , 

for  A  ,4,  .  A*.  ,  4,  he  finds  the  dependence  upon  the  first 

few  f  S  .  One  must  refer  to  the  original  for  details,  especi- 

h 

ally  for  the  scheme  for  approximate  solution  for  the  -f,  i. 


When  tx~0  the  above  analysis  is  precisely  that  for  the 
highest  wave.  The  numerical  results  of  Havelock's  computations 
for  this  case  were  given  in  the  last  section.  He  also  computes 
(fiftl  /  /,  ,  4  (also  /3  for  the  first)  for  two  further 
cases:  e.1"**''  «  O'.  75  and  0.3.  The  agreement  with  results  com¬ 

puted  by  other  methods,  either  those  of  section  27<x  or  similar 
ones,  Is  very  close.  However,  to  establish  the  validity  of  the 
method,  one  must  prove  convergence  of  the  series  . 

The  relation  of  this  method  of  approximation  to  Stokes’ 

"second  method"  (see  section  32(f)  :1s  also  clarified  by  Havelock. 

I 
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f  -  HS~ 


For  this  we  refer  to  the  original  paper. 

3 4 .  Exp lici t  solutions 

Although  it  is  not  in  general  possible  to  give  an  explicit 
exact,  solution  to  a  particular  problem  of  interest,  it  is  possi¬ 
ble  to  give  various  classes  of  exact  solutions  and  then  to  deter¬ 
mine  the  associated  solid  boundaries.  This  is  sometimes  referred 
to  as  an  "inverse  method".  Several  such  methods  for  constructing 
exact  solutions  will  be  discussed  in  section  34  ex.  In  addition, 
there  is  one  periodic  wave  in  infinitely  deep  fluid  which  satis¬ 
fies  the  boundary  conditions  exactly,  the  Gerstner  wave.  This 
will  be  discussed  in  section  34  pi  .  In  section  34f  we  shall  dis¬ 
cuss  briefly  what  may  be  called  pseudo-exact  solutions  due  to 

Davies  and  . . .  In  these  the  e  .act  boundary  condition  is 

replaced  by  a  closely  related  one  which  allows  exact  solution. 
They  derive  their  interest  from  the  fact  that  they  contain  in 
one  family  waves  ranging  from  the  smallest  ampl i tude- length  ratio 
up  to  a  counterpart  of  the  Michel  1  wave.  Furthermore,  the  pro¬ 
cedure  also  can  b©  used  for  pseudo  solitary  and  cnoldal  waves. 
Section  34 2  will  be  devoted  to  an  exact  solution  for  pure  capil¬ 
lary  waves  recently  discovered  by  C rapper  [ 1957 J . 

3 4  ix.  Inverse  methods 

“*'  -  ...  MU  I  T  ■  --■!  II 'll  I  I  1 1  '  -  -----  4 

S  a u treaux ’ s  me t hod .  Possibly  the  earliest  method  capable 
of  generating  a  wide  class  of  steady  irrotational  solutions  is 
due  to  C.  Sautreaux  [1893,  1894,  1901].  It  has  been  rediscovered 
several  times  subsequently,  e.g,,  by  Blaslus  [1910],  Wilton 
[1913],  Richardson  [1920]  and  Lewy  [1952].  F.  Almond  [1929]  has 
given  a  very  comprehensive  treatment  of  the  method  and  of  various 


related  ones.  The  method  may  be  easily  generalised  to  include 
an,  arbitrary  impressed  pressure  distribution  on  the  free  surface 
(see  the  papers  of  Richardson  or  Almond). 

Let  j[  =  X  +  oy,  »  f  ~  cf>*  4  ,  and  take  t  as  the  independent 

variable  The  free  surface  will  be  represented  by  ll//=o  ,  we 
further  assume  <£.  >1  >0  .  In  the  constant-pressure  condition 
on  the  surface,  -jj,  +  J.  ^  -  const.,  it  will  be  convenient  to 

take  the  position  of  the  x-axis  so  that  the  constant  is  Eero, 
and  hence  "J  S  O  .  This  condition  may  then  be  expressed  in  terns 
of  J.ff)  as  follows  [cf  .  (32.50)]: 

Ly?)-T7i]  = 

(34. 1) 

Def lne 

/-rf)  =  r  i’f«-  Jin j 

(34.  2) 

Then  -/*(' f)  =  ,  the  y-coordinate  of  the  free  surface. 

Hence,  from  (34,1) 


j. 

2  = 


yw  yw 

From  (34.2)  and  (34.3)  one  may  now  derive 


-  4/x'V^i  =  Lj'w  +  ywJ 


Elimination  of  J/  b<?tw66n  (34*2)  and  (34<  4)  y io Ids 
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(34.3) 


(34.4) 


(34.5) 


(3d.  6) 


But  then,  since  1/  is  un  analytic  function  of  4  t  at  least 

u 


near 


l/'  »  o 


1 
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'{*)  =  ->'m  +  /a  -  r-: 


<34.7) 


and 


(34,8) 


One  may  now  reverse  the  procedure,  select  an  arbitrary  analytic 
function  p.(4)  satisfying  (34.6)  and  construct  the  function  i  ( f J 
by  means  of  (34.8).  The  resulting  function  will  describe  a  flow 
for  which  ^  ( 9)  Is  the  free  surface.  If  (34.6)  is  satisfied 
only  for  some  range  of  -f  ,  then  for  the  remaining  range  one 
must  treat  the  streamline  s  O  as  a  solid  boundary . 

One  can  use  the  preceding  result  to  construct  a  flow  if  the 
form  of  the  free  surface  is  given.  Let  the  surface  be  given  in 
the  form  *  =  in  a  neighborhood  of  some  point  of  the  surface 

Since  ~  on  the  surface,  we  may  define  jf  ^ h  1 

-  %'($)  ;  CT  is  an  analytic  function  of  /-v  for  real  |(r*.  alt" 

follows  from  the  theorem  of  Lewy  and  Gerber  cited  near  the  begin¬ 
ning  of  section  32  /.  Hence,  from  (34.7), 

rM  *  -tbjrf'-r'1.'  04.9) 


Solving  for  I  f  f*'  ,  one  finds 


fr  -  V^TT77^ 


(34. 10) 


Since  ^  is  also  an  analytic  function  of  ,  the  same  relation 
holds  for  jf  /it*  when  j*.  in  complex,  and  consequently 


It  follows  similarly  from  (34.8)  and  (34.9),  first  for  real  to  , 

/ 

then  for  complex  that 


T  -  *  m  —  ,  C“  j'm)  c  ■ 

*  '  J  1  (34.12) 

Equations  (34.11)  and  (34.12)  thus  provide  a  relation  between  f 
and  determined  by  the  form  of  c r(p)  for  real  a*  . 

Rudzki's  method.  Rudzki  [1898]  has  given  a  different  formula 
for  deriving  exact  solutions.  The  derivation  and  statement  of 
the  formula  below  differ  somewhat  from  Rudzki 's,  but  the  result; 
is  equivalent . 


Let  J  -  iff)  and  write 


/  I  <  9 

Z  s.  -r  e  1=5 

1  >  t  f 


5 


(34. 13) 


The  free -surface  condition  may  be  expressed  as  follows,  from 
(32.61)  and  the  equation  preceding  it, 


* 


Hence 


9  W  1/^  =  0 


'/J 

Sm  $  f  0)  d  ijtj  -for  ^sO( 


(34.  1,4) 


(34.15) 


where  the  branch  <?f  the  cube  root  :1s  taken  which  is  real  for  real 
numbers ,  Combining  (34.15)  with  (34.13)  gives 


}'(*) 


4  @(<p,a) 


f  -  (  fh  9( 4  0)i$] 


This  relation  must  then  . . ild  also  lor  ty'+O  ,  i,e, 


r  , 

£  [  •  im  (9'f  6)  d-f. 


(34. 17) 

As  in  Sautreaux's  method,  we  may  now  reverse  the  above  procedure, 

take  9  1  as  a .  arbitrary  analytic  function  of  i  such  that  6 

is  real  for  +  real,  and  construct  *  VfJ  from  (34,17), 


M-?-1  -JLJlft .thoj.-  *?roni  (34.17)  one  can  derive  Immediately 


a  formula  due 

to  Richardson 

[1920]  for  constructing  exaci 

L  solu- 

tions.  Let 

$'(f)  -  " Sir  S(f) . 

Then  A- V1  - ,  G '  and  ( 

134. 17) 

becomes 

V"(f|  - 
„> 

,  C'tl  r 

j!  0  bi  fr 

;  £  (fj  ^  A- s'1-]. 

(34, 18) 


Again,  inversely,  if  6(f)  is  any  analytic  function  such  that,  for 

:real  *  *  6  .Vf-G'1'1  and  6'  are  real,  (34.  18)  gives  a  corres¬ 

ponding  exact  free -surface  flow. 


:  The  Largest  collections  of  specific  flows  con¬ 
structed  b  i  one  of  the  preceding  methods  are  in  the  paper  of 
Richardson  [1920]  and  a  report  of  Vitousek  [1954],  Several  ex¬ 
amples  are  given  below, 

(1)  In  (34.17)  let  5(f)-  const.  -  a  <0  ,  Further,  take 
the  constant  of  integration  as  zero  even  though  this  results  in 
a  singularity  in  j/  on  the  surface.  One  finds  easily 


$  f  /-Tjsm7  . 

The  free  surface  will  consist  of  only  the 
was  irfi  .  However,  the  ray 


(34. 19) 

ray  unless 

is  also  a  streamline, 


p . ro 


but  not  one  along  which  the  pressure  is  constant  unless  «,  .* - n/6 , 
Hence  it  nuet  fc  i  taken  as  a  solid  boundary  in  general.  The  special 
ease  ot  * ~ is  just  the  flow  (33. 4)  considered  earlier  and  has 
a  corner.  Ono  may,  of  course,  take  any  other  streaml ine 
as  another  solid  boundary  representing  a  bottom,  The  pressure 
remains  positive  everywhere  only  if  -»/<  <  ®i  <  O  .  This  special 
family  of  flows  was  discussed  by  Welngarten  {  1904  ] . 

(2)  If  in  (34.8)  one  takes  j*  {{]*$/&  or  in  (34.18)  takes 
5 iff) •  j /l$fc  4  >  where  c  is  some  fixed  velocity,  one  finds 


£.y(f)  =  -a  -4 


/(^wr 


i 


(34, 


This  yields  a  flow  of  the  sort  shown  in  Figure  48c,  taken  from 
Richardson.  The  infernal  solid  boundary  represents  some  stream¬ 
line  y  z  1^',  < O  „  The  free  surface  corresponds  to  the  segment 

'if''  -  o  ,  O  <  <  C*/l3*  in  the  f-pinne . 

(3)  Let  C  be  some  fixed  velocity  and  let 

G  t)  ZZ  u  3  *  *  _£.  f  5  >  J  «*<  j 

>  1,3  J  > 


in  (34.18).  Then 


«jY«  =  [s+t*.u £(] 

(  Jc  J 

(34.21) 

Here  'f'  =0  corresponds  to  the  free  surface.  The  choice  of  the 
bottom  streamline  is  restricted  by  the  necessity  of  avoiding 
having  the  singularity  at  3  -  tanhj-nt^f  within  the  fluid. 
Figure  48a,  also  from  Richardson,  shows  a  flow  computed  from 
(34.21)  for  5=2  and  csj  . 


I" 
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(4)  Flow  over  a  corrugated  bottom  has  been  Investigate  by 
both  Richardson  and  Rudzki  by  essentially  the  same  method. 

.Pol  lowing1  Richardson,  we  let 


GW  - 


]  a 


i  =»  <  l 


Then 


c\'(  f) 

a 


•  •/ 

/ i 


*ioj  a 


-.,5, ..A;  + 


iC’ 


(34.22) 


Figure  46b  shows  a  flow  computed  from  this  formula  for  3*2,0,; 

(5)  Plows  similar  to  flows  over  a  weir,  under  a  sluice 
through  an  opening,  etc.  have  been  considered  by  a  number  of  the 


cited  authors.  Sautreaux  [1901]  applied  his  formula  (34.8)  with 

-  i  C  »  /,  I 

/  I.  ,,  ,  If  '  * 

M  =  j  C  /i  <i  |  e 


to  obtain  a  number  of  different  flows  of  this 

nature.  Figure  48d  shows  one  of  them.  Lauck  [ 1925]  has  also 
constructed  such  flows.  Richardson  obtained  a  flow  through  an 
opening  by  selecting 


6  it)  = 


3c1 


2 


£  -  &  : 


%  :/l 


Possibly  the  simplest  such  flow,  studied  by  both  Blasius  and 
Vitousek,  is  obtained  by  taking  ^  -  Jc  f/d  in  (34.8);  this  yields 

-J/l 


£  *---•/?  +  f  A/¥  -0 


(34.23) 


The  flow  is  shown  in  Figure  48e. 


Fritz . John  * s . method .  Fritz  John  [  1953 J  has  devised  a  method 


for-  constructing  exact  irrotational  two-dimensional  free-surface 
flows  which  may  be  time-dependent.  Let  denote 
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the  . . . .  velocity  potential.  The  flow  of  particles  cm.  the 

free  surface,  ^  .  will  also  be  described  In  a  Lag  rang  lan 


system 


J  -  ff'K-Oj 


(34.24) 


where  «l  Is  a  real  number  associated  with  a  particular  particle 


F  '( *  -  *ni'  tJj  t ) 


(34.25) 


where  f'  denotes  the  partial  derivative  with  respect  to  y  . 
The  equations  of  motion  (2.7),  reduced  to  two  dimensions  and  to 
motion  along  the  free  surface,  given 


j  y  &  y 
l  fl  «'<«* 


<>  i<-  s  »• 

S1  f1  t  1  ol 


(34.26) 


Since  p  »  const,  on  the  surface,  d -  o  and  (34. 26)  states 
that  hj  Aif1  *t  *  cj,  is  perpendicular  to  e^/3 «  ,  or  that 


{P  <-4*  A  £ 

tt  <r 


r:V.  r)  j 


(34.  27: 


where  rfoytj  is  a,  real  function.  Thus  t  t)  must  satisfy 
the  parabolic  partial  differential  equation  (34.27). 

If  e.  («  j(-)  is  a  solution  of  (34.27)  for  some  r  (<*,  f)  which 


is  real,  for  real  «.  and  if  e  and  £t  are  analytic  functions  of  «• 
and  real  for  real  a.  ,  then  one  may  construct  the  velocity  poten¬ 
tial.  f  i ($,t)  for  a  f roe-boundary  flow  as  follows.  Actually,  we 
shall  construct  F  as  a  function  of  ot  and  t  ,  i.e.  we  shall 
construct  a  function  £  related  to  F  by  ti  ■  For 

real  oc  it  follows  from  (34.25)  that 
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fr  -  S"4~ 


tk 

"ien 


c  f 


m,  f  j 


[o i  fj 


(34.28) 


One  may  now  use  the  last  expression  in  (34.28)  to  extend  analyti¬ 
cally  G  ,  and  hence  G  from  real  to  complex  t  .  By  In¬ 
verting  -s  -&(<*, tj,  one  may  now  construct  (  invertibl  1  i  ty 

follows  from  (2.4)  which  implies  £„.<?«-!  )■ 

It  is  possible  to  prescribe  (*,  t)  and  then  construct  the 
associated  function  j  .  For  it  follows  from  (34.26)  with 

^  ■=  4  fx  (■*,*■),  t)  that 

(34. 29) 

Any  set  of  solutions  >  ^  t)  depending  upon  a  parameter  <*.  yields 
a  function  €. («. t)  defined  by 


e  [>{*  ti,  tj . 


(34.30) 


The  function  r{tn,i)  for  real  «.  is  given  by 


ett  +  t  %■ 


A  €L 


(34.31) 


where  (34.29)  has  been  used  in  obtaining  the  last  expression. 

We  shall  suppose  now  that  the  motion  is  steady  and  make  the 
following  special  choice  of  Lagrangian  parameter  oc  .  Select  some 
fixed  point  J.  0  of  the  surface  (x)  and  for  any  particle  on 

the  surface  let  -  ot  be  the  time  at  which  the  particle  was  at  J.0  . 
Since  the  motion  is  steady,  all  particles  take  the  same  amount 
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of  time  to  travel  from  to  any  give a  point  ^  tad  hence 


fi  («,f) 


e  f 0, ft  tA)  m  £  (t- **) 


1 1  t he n  t o 1 1 ons  f r on  ( 34.27)  t h a t  a Is o 


(34.32) 


r  (  o\  ,t)  s.  v  . 


(34.33) 


Hence  (34.22)  becomes  an  ordinary  differential  equation  in  a 
single  variable,  say  r  s  t  4*  : 


6  #  f  T-)  -  *  r(t)  tfVrJ  *4‘4| 


o 


(34.34) 


It  follows  next  from.  (34.  28)  that  S t) 


£%(a  +  t)  and  thus,  if 


6  f  r )  is  an  anal  y  t  i  c  so  1  u  t  i  on  of  (34 . 3  4 ) ,  re  a  1  f  o  r  re  a  1 


fy  j 

\i  if’ 


IMP* 


p  /  f  y 

'm*<  l«< 


C  !  T 

■mmi>  I  ,|l 


( 3  4 . 3  5  ) 


In  this  case  each  choice  of  a  function  r [ r)  ,  real  for  real 


L  $ 


results  in  a  function  efrj  as  a  solution  of  (34.34),  and  then 
in  a  function  6» (r)  obtained  by  a  quadrature  of  (34.35).  One  may 
invert  jsefc+t)  and  find  F  as  a  function  of  ^  as  in  the  last 
paragraph  or  else  regard 


*1 


<r 


e  ( r) 


jlwx 


r  />. 

a  [  T ) 


Slat**1 


(34.36) 


as  parametric  equations  with  complex  parameter  T  . 

Several  examples  are  considered  by  John, two  of  which  are 


time -dependent 


A  simple  and  interesting  steady  flow  is  obtained 


by  taking  K(x)  =  V  ,  a  constant,  in  (34.34).  Then  (34.34)  and 
(34.35),  after  setting  the  constants  of  integration  equal  to  zero, 


yield 
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■A&vr  F 

V  ) 


*+*  kV!  It  —  z  ^ A 

V **  j  V 


CoS  Vf 


(34.38) 
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The  fra®  surface,  obtained  by  taking  r  real  in  the  first  formu¬ 
la,  Is  a  trochoid#!  curve  without  self-intersections  if  A  <  i/v1 
the  wave  length  is  A  and  the  amplitude  is  A  .  If 

A  <  £/vv»  then  |dF/d}j>0  and  .  However,  JF/.^ 

can  become  infinite  if  d^/dr  -  O  .  Such  points  occur  at 


a 

d 


-L 

i  / 


(34.30) 


In  order  to  avoid  having  them  within  the  fluid,  the  bottom  must 
be  chosen  as  a  streamline  which  passes  above  or  through  these 
points.  Figure  49,  taken  from  John's  paper,  shows  several  profiles 
and  the  associated  streamlines  through  the  branch  points  (34.39) 
computed  for  various  values  ofthe  constant h  when  ;  -  m (this  is 
equivalent  to  graphing  eT / /  for  various  values  of  I  r  /  /,<:  ), 

The  surface  profile  and  bottom  come  closer  together  as  /-i  -~®  *  and 
draw  further  apart  as  0  .  For  A  ■»  0.9  they  are  so  close  that 
they  cannot  be  conveniently  separated  in  the  figure;  in  such 
cases  one  may,  of  course 


FIGURE  49 
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have  reservations  about  the  applicability  of  the  perfect- f luid 
mode  1 . 

The  surface  profile  in  this  example  is  exactly  the  same  as 
in  the  Gerstner  wave  treated  in  the  next  section.  However,  the 
Gerstner  wr ve  is  defined  for  infinite  depth  and  is  not  Irrota- 
tional.  The  flow  described  above  may  also  be  obtained  by 
Sautreaux's  method.  Vitousek  [1954]  has  studied  It  by  this  pro¬ 
cedure  . 

Methods  of  Villat  and  Poncin.  At  the  end  of  section  32  / 
brief  mention  was  made  of  a  pair  of  integral  equations  derived 
by  Villat  [1915]  for  the  determination  of  flows  over  some  given 
bottom  profile  and  with  the  top  profile  also  given  upstream  of 
some  point.  The  method  seems  tu  be  chiefly  useful  as  an  inverse 
method  in  which  the  free  surface  is  giv»n  and  the  other  profiles 
sought.  Villat  has  worked  out  one  case,  but  not  in  complete 
detail,  where  the  top  cover  is  missing  and  the  bottom  has  a  de¬ 
clivity. 

Poncin  [1932]  has  further  elaborated  Villat'w  method  in  the 
direction  of  starting  with  the  fixed  profile  and  finding  the  free 
profile.  Actually,  he  does  not  really  achieve  this.  Instead, 
he  is  able  to  construct  a  flow  for  a  fixed  profile  of  the  same 
general  behavior  as  the  given  one,  but  not  identical  with  it. 

The  method  I1:  applied  to  a  number  of  interesting  special  cases, 
including  flow  over  wavy  bottoms  and  over  bottoms  with  a  decli¬ 
vity.  The  solutions  are  generally  for  large  values  of  the  velo¬ 
city.  The  method  and  results  do  not  lend  themselves  to  a  brief 


summary . 
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34 ft  ■  Gerstner ' s  wave 

Gerstner’s  wave  [1802]  is  apparently  the  first  flow  to  have 
been  discovered  which  satisfies  exactly  the  condition  of  constant 
pressure  on  the  surface,  and  is,  In  fact,  one  of  the  earliest 
papers  on  water-wave  theory.  It  was  subsequently  rediscovered 
by  Rankine  [1863],  As  will  be  shown  below,  the  motion  is  not 
i rrotational .  This  fact  itself  would  not  be  enough  to  rule  it 
out  as  a  mathematical  model  for  real  periodic  waves.  However, 
the  direction  of  the  vorticity  is  such  that  it  is  difficult  to 
conceive  of  a  scheme  whereby  such  a  wave  could  be  generated  in 
nature. 

The  motion  is  most  easily  described  in  Lagrangian  coordi¬ 
nates.  Each  particle  is  associated  with  a  pair  of  parameters 
(a,b),  b  iO  .  However,  (a,b)  does  not  represent  a  particular 
position  of  the  particle  at  some  time  t0  ,  but  instead  a  mean 
position.  Hence,  instead  of  (2.3)  and  (2.4)  we  need  require  in¬ 
steady  only  that,  the  determinant  Jj  of  those  formulas  be  inde¬ 
pendent  of  t ,  The  motion  is  described  by  the  equations 


If  £- C  is  taken  as  the  free  surface,  the  motion  evidently  re¬ 
presents  a  wave  moving  to  the  left  with  velocity  o  =  c r/Wi  ,  while 
the  particles  themselves  describe  in  a  counter-clockwise  direc¬ 


tion  circular  paths  about  the  points  (a,b)  associated  with  the 
particles.  The  surface  A  -  O  is  a  trochoid  and,  in  fact,  each 
curve  Am  const.  <0  is  also  a  trochoid.  In  order  that  there 


shall  be  no  self-intersections,  one  must  have 
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A  '  •  (3  4.41) 

lit  order  to  verify  that  the  motion  is  kinematically  possible,  it 
is  necessary  to  show,  as  noted  above,  only  that  the  Jacobian 
'•%,)/&(€>  4)  Is  independent  of  t.  An  easy  computation  shows 
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(34.42) 


so  that  the  continuity  condition  is  satisfied.  Next  one  must 
show  that  the  pressure  is  constant  along  the  free  surface.  We 
shall,  In  fact,  show  more,  namely,  that  it  is  constant  along  any 

line  A  -  const,  c  0  ,  provided  .  To  see  this,  introduce 

the  equation  (34.40)  Into  the  first  of  equations  (2.7)  A 
straightforward  computation  yields 


A 


(  Z  T~ 


Or  Cr 


a  i  r  { iv 


J  kVi  CX  •+  0“  t  J 
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If  the  pressure  Is  constant  along  the  surface,  then  ®b/6(X-Q, 
This  can  only  hold  if 


crv  =  ^  ^  (34.44) 

However,  if  c*v  -  <j.  **  ,  then  cl  p  /  ia.  for  all  4r  ,  so  that  each 

curve  $r m  const,  is  an  isobaric  curve.  Although  we  shall  veri  fy 
this  fact  directly,  it  now  follows  immediately  from  Burnside’s 
theorem  in  section  32  fi  that  the  motion  cannot  be  irrotational , 

A  direct  computation  of  the  vorticity  vector  is  facilitated  by 
noting  that 

u  •=*  ■*-  =  A  re~'^cos  +  rt)  =  - 

/\cr  ^  d  +  erf)  =  <T  (  >  -  o)  , 

Zt 


Vr 


(34.45) 


The  right-hand  side  of  (34.46''  may  be  computed  from  (34  40)  by 
application  of  the  rules  of  inversion  for  partial  derivatives. 
The  f ina 1  resul t  is 


~d  v 
>  x 


x. . 

rV  h 


*  A 


,  z~x- 


£ 


(34.47) 


the  negative  sign  indicates  that  the  vorticity  is  directed  oppo¬ 
sitely  to  the  orbital  notion  of  the  particles. 

We  shall  omit  a  discussion  of  the  construction  of  the  curves 
Jfj~  .  const.,  streamlines  in  a  coordinate  system  moving  with  the 
waves;  it  may  be  found  in  Lamb  [1932,  S  251],  Mi lne -Thomson  [195(3, 

(5  14.81]  and  in  Kochin,  Kibel’  and  Roze  [  1948,  ch.  8,  ^  16]  to¬ 
gether  with  reproductions  of  Gerstner’s  original  curves.  It  is, 
however,  of  interest  to  note  that  there  is,  according  to  (34.41), 
a  "highest"  wave  dL  ratio  lA/X  ~  f  ,/ff  -  0.318,  a  figure  which  may  be 
compared  with  the  value  0.142  for  Michell's  wave.  The  highest 
Gerstner  wave  has  a  cusp  at  the  crests,  a  further  indication  that 
the  motion  cannot  be  < rrotational . 

The  pressure  distribution  may  be  found  by  substituting  (34,4)) 
in  the  second  equation  of  (2.7),  using  (34.44),  and  integrating. 

The  result  is 


r  n  "  a 


/  1  '*  & 

’ 1  p 


~  f>  cr  *  A  if 

(34.48) 

A  computation  of  the  potential  and  kinetic  energies  over  a  wave 
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length  yields 


T  =  V  =  ^  A~[  l  - 


(34. 49) 


Finally,  we  note  that  nowhere  in  the  preceding  analysis  have 
we  made  use  of  homogeneity  of  the  fluid,  i.e,  the  Gerstner  wave 
also  represents  an  exact  solution  for  an  arbitrary  heterogeneous 
fluid  (with  constant  along  streamlines).  Moreover,,  Dubreil- 
Jacot in  [1935]  has  shown  that  the  Gerstner  wave  is  unique  In  this 
respect . 

Gerstner 's  wave  is  defined  only  for  infinite  depth.  One 
may  ask  if  a  similar  wave  exists  for  finite  depth.  "SlW  lar" , 
in  this  context,  will  be  taken  to  mean  a  periodic  wave  which 
satisfies  exactly  the  constant-pressure  condition  on  the  free 
surface  and  for  which  the  particle  orbits  are  closed.  Dubreil- 
Jacotin  [1934]  has  proved  the  existence  of  such  a  wave  and  showed 
that  it  is  unique  when  the  period  is  fixed.  However,  this  motion 
cannot  be  given  explicitly  except  in  the  case  of  infinite  depth. 
Methods  of  approximate  computation  of  the  wave  have  been  given 
by  Kravtchenko  and  Daubert  [1957]  and  Gouyon  [1958]. 

3 4  jf'  .  Pscudo-exact  solutions 

Although  the  solutions  of  this  section  are  not  really  solu¬ 
tions  satisfying  the  exact  boundary  conditions  formulated  earlier, 
they  are  exact  solutions  to  a  closely  related  problem,  also  with 
a  nonlinear  boundary  condition.  Their  interest  derives  from  the 
fact  that  it  is  possible  to  encompass  within'  one  explicit  for¬ 
mula  waves  of  all  amplitudes  up  to  a  highest  wave  analogous  to 
Michell’s  wave.  The  procedure  also  allow®  explicit  construction 
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of  solitary  and  cnoidal  waves.  It  is  possibly  a  misnomer  to 
call  these  solutions  pseudo-exact,  for  one  may  also  interpret 
them  as  the  first  term  in  a  certain  series  solution  of  the  cor¬ 
rectly  posed  problem.  In  this  sense  they  are  analogous  to 
Havelock's  approximation  procedure  described  in  section  23/2  . 

The  work  to  be  described  has  appeared  in  a  series  of  papers  by 
Davies  [1951,  1952],  Packhan.  [1952]  and  Goody  and  Pnvies  [  1957], 
The  motion  will  be  described  in  terms  of  the  variables  In¬ 
troduced  in  (32. 85) ,  CJ  —  0  +  The  alteration  in  the  ooundary 

condition  consists  in  replacing  (32.89)  by 

^Y  r-  (34.50) 

where  /  is  some  fixed  constant  chosen  so  that  niB  is  a  good 
approximation  to  Z  f  '?  If  one  wishes  to  consider  (34.50)  as 

the  first  term  in  a  series  approximation  to  (32.89),  one  may 
expand  s>r  S  in  a  series  in  c.  *■  and  express  (32.89)  as 


ie 


y .  3  5  +• 


_4 

S\ 


3 

J  i  ►»,  )  &  -+ 


(34.51) 


In  this  case  (34,50)  with  A-  Vi  represents  the  approximation 
obtained  by  keeping  only  the  first  term  of  (34.51).  However, 
we  shall  not  pursue  the  approximation  procedure  and  refer  to 
Davies  [1951]  for  further  information.  It  will  be  convenient 
to  reformulate  the  boundary  condition  (34,50)  as  follows: 


in 


d.  -f ' 


£  /  -  0  'for  f-  0 


(34.52) 
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*<3m> 


or,  after  introducing  the  new  variable  >((f j  s  6  "  s uz/c1*  as 

Ivy,  -4  +  i  =  o  +**  o , 

I  \  **■  CJ  /J 


<34 . 53 ) 


In  order  to  proceed  further,  we  must  further  specify  the 


nature  of  the  wave  motion.  Let  us  suppose  the  notion  to  be  peri¬ 
odic  with  wave  length  X  and  take  the  fluid  to  be  Infinitely  deep. 
We  again  introduce  the  f-plane  of  (32.91)  and  take  coordinates 
as  in  Figure  46.  The  expression  in  curly  brackets  in  (34.53)  is 
a  regular  analytic  function  of  f  inside  the  unit  disc  of  the 
I  ~P Lane  with  vanishing  imaginary  part  on  the  boundary,  hence 
is  a  constant.  Since,  for  £  =  O  (  i.e.  as  "f"  --»-«=») ,  >'  =  <  and  -r, 

the  constant  must  be  3i«/cJ  .  Thus  \  must  satisfy  the  differen¬ 
tial  equation  - 


ire  cJ 


(34,54) 


The  solution  is  easily  found  to  be 

-,3  Ut/c3 

^  =.  |  -+  A  e 


(34.55) 


Referring  to  Figure  46,  we  aee  that,  if  {  =  0  ,  ^  r  %*  / c  >  where 


is  the  absolute  velocity  at  a  crest.  Hence 


3 


A  *  -%  -  / 

c 


(34.56) 


Sine©  9  must  also  vanish  at  {pzi jcX,  i.e.  the  left-hand  side 
of  (34.55)  must  be  real,  we  must  also  have  (3$ jA^Jj <■  ^  s  T  ,  or 


v 


3  I'  ^  / 1 7! 


(34.57) 
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Mot©  that  if  )z  I/s  ,  the  relation  between  c/  and  a  is  the 


same  as  in  the  infinitesimal-wave  theory.  The  solution  (34.55) 


may  now  be  put  into  the  following  form: 


=  i-(<- 1= 


>  .•  f  / 

-ii-n+rt, i 


(34.58) 


where  0  <  5  „ < 


If  ire  ,  then  ur  =  c  and  the  flow  is  uniform. 


O  ,  then 


j  f  -  J  /c  )  i 

=  e  I  1  -  e 


(34. 59) 


and  near  I  -  0,  t  £.  A  ,  ii'-A  there  is  a  corner  in  the  wave 

profile  with  the  two  tangents  making  the  same  angle  120°  as  in 
Stokes'  theorem  (near  -f  -0  ,  (33.5)  gives  ur*  r  <  j-  9  |  ,  (34.58) 

gives  ur3=-*3i^l  ).  Hence  this  wave  corresponds  to  Michel! ’s 


highest  periodic  wave.  The  ratio  of  amplitude  to  length  of  this 


wave  may  be  computed  from  the  following  expression  for  the 


trough: 


)k  -  A  ft. 


<-  .  L 

0 


rc)  f. 

I  -‘i*#/ S, 
I  J  —  *  M 


By  expanding  in  a  series  and  integrating  term  by  term,  one  finds 


—  -  O.  Ill, 

> 


(34.60) 


We  recall  that  the  value  for  Michell's  wave  was  0.142. 


If  the  depth  of  fluid  is  finite,  one  must  add  the  additional 


boundary  condition  . Mx!  toT-^z-o  ,  as  well  as  for  fl-0 
and  i  5  if  the  motion  is  to  be  periodic.  The  determination 


of  )(  now  becomes  too  difficult  to  carry  through  briefly.  How- 
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ever,  an  explicit  solution  is  still  possible  and  has  been  worked 
out  by  Davies  [1952]  and  further  investigated  by  Goody  and  Davies 
[1957],  Similarly,  a  "solitary  wave”  can  be  explicitly  constructed 
which  satisfies  the  boundary  conditions  1~  ' /\  -  O  for  l/'  -  -  fJ  and 
for  -  o  ,  o  >  and  '/  )  as  $  -*  t  .  This  has  been  done 

by  Pack ham  [1952],  Either  of  these  problems  leads  to  the  follow¬ 
ing  differential-difference  equation  foi  )(•  *)  : 


^  (  f  ■+  <  o) 


r  * 


yj 


x 


(34.61) 


it  may  be  established  in  a  manner  similar  to  that  used  in  deriv¬ 
ing . (22,30)  or  (32.80). 

34  I  .  Pu re  capillary  wa v e s 

The  first  investigation  of  periodic  progressive  capillary 
waves  satisfying  the  exact  boundary  condition  is  apparently  due 
to  N.  A.  Slezkin  [1937],  He  formulated  the  boundary-value  prob¬ 
lem  In  the  same  manner  as  will  be  done  below,  reduced  it  to  solu¬ 
tion  of  a  nonlinear  integral  equation  analogous  to  Nekrasov's, 
and  proved  existence  and  uniqueness  of  a  solution.  However,  he 
apparently  did  not  observe  that  an  explicit  solution  was  possible 
for  infinite  depth  of  fluid.  This  was  discovered  by  Crapper 
[1957],  following  a  different  .and,  in  fact,  more  elementary  method. 

We  shall  consider  the  motion  as  a  steady  one  in  which  the 
fluid  moves  to  the  right  with  velocity  c  as  ^  -  sc  .  The  exis¬ 
tence  of  a  complex  velocity  potential  f  will  be  assumed 

and  the  free  surface  ^  s  yfx)  will  be  taken  to  correspond  to  the 
streamline  1^s0  as  usual,  It  will  also  be  convenient  to  make 


till . . . . . . 
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u»®  of  the  variable  u>«a?**r  introduced  In  (32.85).  If  h0  is 
atmospheric  pressure,  then  from  Bernoulli's  integral 


l’  i  r  \  =  r.  ♦  i  r 


(34.62) 

(we  recall  that  gravity  is  being  neglected).  The  dynamical  con¬ 
dition  at  the  free  surface  [see  (3.8)  and  (3.9)]  is 


J  P  0  "  “* 

1% 


ff  »  y  'ml 


i/J. 


(34, 63) 

Before  combining  (34.62)  and  (34.63),  we  recall  that  the  curvature 
of  a  streamline  at  any  of  its  points  is  given  by  J8/J4.  where  *&■ 
is  arc  length  along  the  streamline.  Hence,  we  may  combine  (34.62) 
and  (34.63)  to  obtain  the  following  boundary  condition 


jrU'-'.')  =t  -  t if  il  i.r  +  = 

L  “■  a  a-  *  h 


(34.64) 


or 


XI  f  i-  -  % 
11  1  f  c  J 


CJ& 


fOr  ^  —O  . 


(34.65) 


Since  ^.*c-eT  and  since  S <9/7)  f  -  s> r/S  f  from  the  Cauchy-Riemann 
equations,  (34.65)  may  be  written 


~br 

If 


Y=(e~t-e^)  =  -  y  s  U *  Ur  f 


0. 


(34.66) 

The  problem  is  now  to  find  a  function  to(-f)  analytic  for  V f'iD  , 
such  that  as  and  such  that  the  imaginary  part  T 

satisfies  (34.66).  However,  since  the  boundary  condition  (34.66) 
involves  only  T  unlike  its  analogue  (32.89)  for  pure  gravity 
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waves,  it  Is  possible  to  solve  first  for  the  harmonic  function 
and  then  to  find  <9  later. 


We  as  suite  that  a,  solution  can  be  found  'which  satisfies 


"d'T 


e)  ^ 


r 


* 

T 


and  proceed  to  verify  the  assumption. 


(34. 67) 

I  n  t  e  g  rati  n g;  ( 3  4 . 67 )  ,  we 


obtain 


j  fi*if  r 

J.  » 


Uf 


>  j  ^(4>) 


{34. C 8) 


where  H  7V/)  =  IX"/'')  and  G{rpj  is  an  arbitrary  function,  or 


r 


,  /  •>  -  r '  4 1 

**  . .  -n  . _ — - . 


(34. 69) 


Since  T  is  a  harmonic  function  of  <f>  and  ‘/X  ,  Laplace's  equation 
must  be  satisfied  by  (34.69).  This  yields  an  equation  to  be  satis 
fled  by  X  and  ,  in  which  the  two  variables  can  be  separated.  We 
shall  not  repeat  the  detailed  analysis,  which  is  typical  of  that 
occurring  in  separation-of -variables  problems.  The  final  result 
is  that  /  and  /  must  satisfy 


X'L  =-  a,  *■  ’-l  /'  «  - 

Y'v  -‘»1Yv-a3  Y\ 


(34.70) 


where  <3,  ,  av  ,  Qj  are  arbitrary  constants.  Crapper  states  that 
the  full  equations  may  be  used  to  construct  a  solution  for  fluidof 
finite  depth,  but  that  it  is  sufficient  to  set  =  O  for  infinite 
depth  (in  view  of  Slezkin  ’s  result,,  this  is  presumably  also  neces¬ 
sary).  Since  T  is  real,  we  shall  also  take  X  and  T  to  be  real. 
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If  one  does  set  ftjSS  and  assumes  <0  ,  at  >o  ,  the  following 
give  real  solutions  of  (34.70) 


xfy)-/'-i  ,m=f' X 


/  > 


(34.71) 


where  ,t  and  b  are  real  constants.  A  glance  at  (34.69)  shows 

that  t  is  Independent  of  the  choice  of  0. 

'1 


It  will  be  con¬ 


venient  to  let  -  (v  / c k  ,  where  iv  >  o  .  One  may  determine  £ 
from  (34.07),  for 


tr 


tii'1'  tZ 


or 


i  6  7  /  <»  £  '  ■+  1  -  ... 

Pl  —  . —  -l! _ —  —  "3 

(34.72) 

Since  £•  is  to  be  real,  we  must  evidently  have 

!*I  ^  j 

p  Cx  '  '  (34.73) 

Since  if"  adds  only  a  real  constant  to  p  we  may  select  it  at  our 
convenience;  we  take  F  =  0  .  Substitution  of  (34.71)  into  (34.69) 
gives 

_  -j  ,  coS.»(*v  Wc.  +£  )  +  tos 

T  —  J  O  Q  - ; - - - 

J  C  05 (  wV'/  c  +  £  J  -coS'V  <^/c  j 


c  os  (i  f'  /Vc  *+(£’)  +  CoS (*+.  ii/c  J 

. — !**■■-- - - - - —  i  I,  .  .  _  _  __  

Col(7*»yVt  <M  £  J  -  C  35 


=■  !<*  [cot  +;£jtot 


(34.74) 
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The  analytic  function  m>  ,  which  has  v  as  imaginary  part  and 
wnich  approaches  aero  as  iL/'-%~ao  ,  is  given  by 


le  then  have 


(34,75) 


—  -  c  fit  r  -  w  c  st  lltrt/c  ■**£]  =  c  r  lf  **  Vc  -fj 

I  .  1  V.  '  * 


(34. 76) 


From  this  one  may  solve  for  j,  in  terms  of  f 


f 


a  T 


f 


2  c 

VV\ 

ta"  l(  t"  >EJ 

•+  C  0  A  it 

....  <  ±L 

1 

-+• 

-  f  *  ^  (/c.  - 1  J 

•+*  £  '■ 

-  4 

/ 

-+  ;  1b 

'  +  3e 

(34. 77) 


where  the  constant  has  been  chosen  so  as  to  make  c ^  reduce  to  t 
when  B  =  O  .  It  i  evident  that 


"iJLy  ^ 

*<Vs 


J  =  }(*) 


21 


so  that  the  streamlines  are  periodic  in  the  x-direction  with  wave 
length  A  =  LtT/U*  , 

The  surface  streamline  is  obtained  by  setting  1 ^~o  .  After 
separation  of  real  and  Imaginary  parts  in  (34.77)  the  equation  for 
the  surface  becomes: 

3  Sha  f*ji/c 
f  -+  B  L  -t  1 2>  c*1  **  'PA 

f  -t  B>  c-os  4/c 
t  *+  B'1'  +  23  CoS  kw  4/c 

(34.78) 
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with  $  serving:  as  a  parameter*  There  is  a  crest  when  fl5>  *  o  and 


a  trough  when  6 


'j  ii  'L  Hl 


.i*1 


The  di flare ace  in  the  values  of  %  yields 


the  following  expression  for  the  ratio  of  total  amplitude  to  wave 


ength : 


edit*!' 


4  3 

t  t  /  —  pjs 


if 


>1 
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(34.79) 


Equation  (34.72)  then  provides  a  relation  between  A  />  and  iwT/f  c 


V 


•which  we  may  write,  for  example,  as 


C 


i 


m 


bVi  / 


X 


if  / 

)  Yu 


p 


f  -4  *v"i 

'  t(m 


/ 


-i—  a, 
£4 


i.  V 


(34.30) 
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If  this  formula  is  compared  with  (27.29),  if  should  be  kept  in 


ind  that  0..  is  here  the  total  amplitude  and  that  in  (27.29)  /, 


is  a  length  associated  with  the  half  amplitude.  The  formulas  are 


cons is tan t 


As  a/),  increases,  the  surface  profile  becomes  steeper  and 


steeper  near  the  troughs  until  the  two  sides  finally  touch.  This 


occurs  for  -  0.730 


A  wave  of  these  proportions  may  be  con 


3ldered  as  a  "highest"  capillary  wave,  an  analogue  of  Michell’s 


wave,  although  the  nature  of  the  limitation  is  different.  Figure 


50,  reproduced  from  Cropper’s  paper,  shows  the  profile  of  thi 


!*. 


wav©  together  with  other  streamlines.  It  is  a  consequence  of  the 


form  of  the  dependence  in  (34.77)  that  the  other  streamlines  in 


Figure  50  may  also  serve  as  surface  profiles  for  different  value 


ra 


Q'f  &*/\  f  i„e«  for  different  values  of  J3 


It  is  not  surprising 


P 


of  course,  that  the  profiles  are  similar  to  the  middle  three  of 


I* 
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3  5 , 


Existence 


t heorems 


In  the  various  applications  of  the  approximate  theories  of 
chapter  D  and  E  It  is  tacitly  assumed  that  there  is  an  exact  solu¬ 
tion  which  is  being  approximated.  Without  knowledge  of  conditions 
for  existence  and  uniqueness  of  a  solution  to  a  particular  prob- 
le,,1»  the  status  of  an  approximate  solution  is  somewhat  anomalous 
and  one  must  rely  upon  comparison  with  experimental  results  for 
conviction  concerning  the  correctness  of  the  solution.  However, 
such  comparison  Is  not  a  satisfactory  criterion,  for  in  the  ori¬ 
ginal  formulation  of  a  problem  one  will  usually  have  already  made 
a  decision  about  the  mathematical  model  of  a  fluid  which  will  be 
used.  Thus,  if  one  has  assumed  a  perfect  fluid  (as  we  usually 
have>  and  . . .  "ad*  »  further  mathematical  approximation  in  solv¬ 

ing  the  problem  at  hand,  the  validity  of  this  approximate  solution 
must  first  be  established  before  comparison  of  the  predicted  re¬ 
sults  with  experimental  measurements  can  be  used  as  a  criterion 
of  the  applicability  of  the  fluid  model.  Without  this  additional 
knowledge,  the  comparison  of  approximate  solutions  with  experi¬ 
mental  results  must  be  considered  in  some  sense  to  be  second  best 


even  though  valuable  evidence  may  be  provided  by  good  agreement 


in  a  wide  variety  of  situations.  . 

Uni ortunate 1 y ,  existence  and  uniqueness  proofs  in  exact 
water-wave  theory  have  generally  been  difficult  to  establish,  and 
have  usually  been  obtained  for  only  rather  restricted,  although 
physically  important,  situations.  Many  of  them  are  very  recent 
and  some  rely  upon  methods  of  functional  or  topological  analysis 
which  cannot  be  briefly  expounded.  Although  some  proofs  are  so 
constructed  that  approximation  methods  are  inherent  in  them, 

. 
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others  are  only  able  to  assert  the  existence  of  a  solution  with 


no  indication  of  how  to  obtain  it  approximately,  Proofs  are 
still  lacking  for  many  relatively  simple  but  important  problems, 


for  example,  Mlchell's  highest  wave  and  standing  water  waves. 

No  attempt  will  be  made  to  give  an  exposition  of  the  mathe¬ 
matical  methods  which  have  been  used  in  establishing  the  various 
existing  theorems.  Instead  only  a  discursive  account  will  be 


given 


the  nature  and  limitations  of  the  known  theorems 


35ot .  Irrotationnl  waves  -  infinite 
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Proof  of  the  existence  of  periodic  waves  of  permanent  type 
in  infinitely  deep  water  was  first  given,  by  Nekrasov  [  1921,  1922) 
i  n  a  j  ou  rn  a  1.  of  ve  r  y  res  trie  ted  d  i  s  t  r  1L  bu  t  il  on ,  Sh  o  r  1 1  y  t;  he  re  a  f  te  r 
Levi-Ci vita  [1925]  gave  another  proof  along'  quite  independent 
lines.  Further  proofs  were  given  by  Neumann  [1929]  and  by 
Lichtenstein  [1931],  these  being  more  closely  related  to  Nekrasov’s. 
A  new  treatment  of  Levi-Ci vita ’ s  proof,  due  to  Littman  and 


Nirenberg,  is  contained  in  Stokers1  Water  waves  [1957,$  12.2], 
.Also,  Nekrasov  [1951]  has  recently  published  his  researches  in  a 

m  o  r  e  a  c  c  e  s  s  i  b  1  e  f  o  rtu . 

,«» 

Nekrasov’s  method  requires  proving  that  there  exists  a  solu¬ 
tion  0(y)  to  his  nonlinear  integral  equation  (32.104).  His  pro¬ 
cedure,  in  brief,  is  to  assume  an  expansion  of  $W)  in  powers 


of  the  parameter  -  t*  -  3  >  0  , 


(35. 1) 


then  to  derive  equations  relating  each  <9.  to  ones  of  lower  index, 

¥ 

and  finally  to  show  that  the  series  converges,  H  -3  ls  chosen 
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dLjl* 


as  a  starting’  point;  because  It  is  the  first  eigenvalue  of  the 
"linearized"  equation  (32.104),  i.e.  the  one  obtained  by  replac¬ 
ing  the  quotient  containing  3  by  simply  m  <9(/i)  .  This  cor¬ 
responds  to  the  infinitesimal-wave  theory.  Proof  of  convergence 

requires  that  m/  be  sufficiently  small,  and  positive,  but  no 

l 

estimate  of  radius  of  convergence  is  obtained.  On  the  other  hand, 
the  method  does  allow  computation  of  explicit  approximate  formulas 
for  quantities  of  interest. 

Levi-Civi ta  also  works  with  the  variable  go  ,  treating  It 
as  a  function  of  the  variable  [  introduced  in  (32.90).  Hence 
his  formulation  of  the  problem  is  essentially  the  same  as 
Nekrasov's,  l.e.  he  is  seeking  a  function  ,  regular  in  the 

disc  I  j’  |  <  i  ,  vanishing  at  y  =0  and  satisfying  (32.97)  on 
If)  =  1  and  some  further  condition  assuring  that  \*r/  c  -  j  ■  <■  \  c  I  . 
His  procedure  for  finding  such  a  function  is  to  expand  both  to 
and  5  (  -  J  A  / 2  tt  C  in  a  power  series  in  a  parameter  (M  O  ; 


go  i- 


«Q 

g—  /  '}  V*i 

y  to  { ( i  h 

i -  h  «  I 

h  s  I 


(35.2) 


where  the  functions  (  Q  and  the  constants  are  to  be  deter¬ 
mined  by  the  boundary  conditions.  The  first  terms,  =  -  j  J  , 

\  ~  O  ,  correspond  to  infinitesimal  waves,  so  that  the  para¬ 
meter  <-*.  is  essentially  the  amplitude /wave ~lengfh  of  this  approxi¬ 
mation.  Levi-Civita  establishes  the  convergence  of  the  series 


iHu, 

(35.2)  for  sufficiently  small  values  of  yH  .  No  estimate  of  a 
radius  of  convergence  is  given,  but  Hunt;  [1953]  has  stated  that 

'din 

an  examination  and  refinement  of  Levi-Civita's  inequalities  show 
that  convergence  is  established  for  amplitude -wave  length  ratios 
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u-3  to  1/98.  The  procedure  lends  Itself  to  explicit  computation 
of  higher- order  computations,  and,  in  fact,  he  carries  them  out 

through  '  -  5.  I . . . Civita  further  derives  the  interesting  theorem 

that  irrotational  waves  of  permanent  type  must  be  symmetric  about 
vertical  lines  through  crest  and  trough.  Nekrasov  assumed  this 
at  the  outset. 

Neumann  and  Lichtenstein  (the  latter's  approach  is  simpler) 
derive  a  coupled  pair  of  nonlinear  integral  equations  and  put  them 
into  a  lorm  such  that  Schmidt's  theory  of  nonlinear  integral  equa¬ 
tions  is  applicable.  Iterative  methods  of  solution  can  be  used 
to  obtain  approximate  formulas. 

J >3 ^  .  I rrotat  ion?  i  waves  -  horizontal  &o]j  (gm 

When  the  fluid  is  infinitely  deep  and  the  notion  periodic, 
the  only  Independent  dimensionless  parameter  besides  the  arnplitude- 
wave  length  ratio  is  r1/  :  >  .  When  the  fluid  is  bounded  below 

b>  a  horizontal  plane  at  mean  depth  r  ,  then  a  new  parameter , 

HI  /  A  I 

say  £  /  £  h  ,  must  enter  into  any  solution.  However,  other  in- 
dept,  t  lent  sets  of  parameters  may  be  used,  and,  in  particular,  dif¬ 
ferent  choices  of  a  perturbation  parameter  have  led  earlier  to 
different  approximate  solutions  for  finite  depth.  Thus,  in  sec¬ 
tions  1 4  fl  and  27  one  finds  the  first  and  higher  approximations 
for  periodic  waves  of  permanent  type  when  A/h  Is  taken  as  a 
perturbation  parameter,  whereas  in  section  31  one  finds  approxi¬ 
mations  to  two  further  types  of  waves  of  permanent  type,  one  of 
them  periodic,  corresponding  to  a  different  choice  of  parameter 
and  a  different  method  of  approximation.  In  each  of  these  cases 
there  arises  the  question  as  to  whether  there  exist  waves  of 


permanent  type  satisfying  the  exact  boundary  conditions  for  which 
these  waves  may  be  considered  approximations .  In  each  case  the 
answer  is  affirmative, 

Waves  of  smal . 3  amplitude.  The  first  proof  of  the  existence 

of  periodic  progressive  waves  in  fluid  of  finite  depth  ts  due  to 
Struik  (  1926 1  His  method  of  analysis  is  similar  to  Lev l -Civi ta  '  s 
for  infinite  depth .  Existence  of  the  desired  wave  is  established 
for  each  value  of  -  /*  «  <  !  and  for  each  sufficiently  small 
value  of  A.  /  /  ,  where  the  bound  on  A  /  A  depends  upon  c  /o,  k 
Hunt  [1953]  has  recently  corrected  some  errors  in  the  proof  which 
d  i  d  n  o  t;  i  n  v  a  1  :ii  d  a  t  e  it  b  u  t;  which  re  s  u  1 1  e  d  i.  n  i  n  c  o  r  r  e  c  t  a  p  p  r  o  x  i  - 
mate  formulas,  11 

Nekrasov  [  1928,  1951]  was  also  able  to  show  that  his  integral, 
equation  for  £  ,  as'  modi  fled  for  finite  depth  [see  (32.104)  and 

(32.106)],  had  a  solution,  thus  providing;  an  independent  proof. 

As  was  the  case  for  infinite  depth,  Nekrasov  assumes  that  the 
waves  are  symmetric  about  verticals  through  trough  and  crest : 

Struik  proves  this.  Krasnose 1 sk i 1  [1956]  has  recently  applied 
topological  methods  of  analysis  to  Nekrasov's  equation  and  es¬ 
tablished  not  only  existence  of  solutions  for  m,  in  the  neighbor¬ 
hood  of  the  eigenvalues  of  the  linearized  equation,  but;  also 
their  uniqueness  and  continuous  dependence  upon  m.  . 

§?A±±*jy  and  cnoldal  waves,  Lavrent’ev  [1943,  1947]  was  the 
first  one  to  establish  the  existence  of  cnoidal  and  solitary 
waves.  Cnoidal  waves  are  not  mentioned  by  him  by  name ,  but,  in 
fact,  their  existence  for  sufficiently  large  wavelength  is  es¬ 
tablished  along  with  that  of  the  solitary  wave,  the  latter  being 


/  / 


obtained  as  a  limiting  case.  The  detailed  exposition  of  the  re¬ 
sults  [1947]  is  unfortunately  both  difficult  of  access  and  diffi¬ 
cult  to  read,  and  relies  upon  earlier  theorems  of  the  authors, 
Although  the  perturbation  parameter  appears  at  first  glance  to 
be  taken  as  r  "  =  -  1  -  i  r  ■  'J  ,  which  for  the  solitary  wave  would 
be  in  contradiction  with  (32.52),  the  quantity  r  is  not  really 
n;ean  depth  but  a  related  quantity  which  varies  with  the  wave 
length  of  the  approximating  periodic  wave.  Friedrichs  and  flyers 
[1954]  by  a  completely  different  procedure  have  established  the 
existence  of  the  solitary  wave.  Their  perturbation  parameter  is 

‘1*  \  t  "...  -  ,  ,  h\  ^ 

essentially  £  *1->1  (actually  it  is  fit*  =  -  j  )  }. 

The  point  of  departure  is  again  the  boundary  condition  (32.89) 
for  the  func t  ton  .  However,  an  Integral  equation  is  formu¬ 

lated,  then  altered  by  a  change  of  variable  cf  =  a  <p  ,  'j/  -  V  .  The 
different  rates  of  stretching  correspond  to  those  of  section  10/3  . 
(Something  like  this  also  occurs  in  Lavrent'ev's  proof,  but  is 
disguised  in  his  theorems  on  conformal  mapping  of  narrow  strip¬ 
like  regions . )  An  iterative  procedure  is  used  to  prove  existence 

% 

of  a  solution  for  sufficiently  small  values  of  S 

Littman  [1957]  has  used  a  method  somewhat  similar  to  that 
of  Friedrichs  and  Hyers  to  establish  the  existence  of  cnoidal 


waves  satisfying  the  exact  boundary  conditions.  However,  as  a 
parameter  he  has  used  essentially  r / }  ,  where  h  is  the  mean 

depth  and  A  the  rave  length.  It  is  demonstrated  that  solutions 


exist  for  values  of  C  Aj.  k  which  are  both  greater  and  less  than 


1.  Figure  51,  modified  slightly  from  one  in  Liftman's  paper 
shows  In  a  qualitative  fashion  the  relation  between  the  dimen* 
sion less  parameters .  The  dotted  lines  enclose  values  of  the 


. . . .  !j|!HttUiiwjit'iu . . 


parameters,  again  in,  a  purely  qualitative  way,  for  which  solu¬ 
tions  have  been  demonstrated  to  exist.  (Here  A  is  the  modulus 


of  the  elliptic  functions  and  K  Is  the  complete  elliptic  ii 


tegral  of  the  first  kind.  <n  serves  as  a  parameter  In  certain 


approximate  for mu 1  as . ) 
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FIGURE  51 


35/  .  I r rotational  waves  -  other  conf igurati 


Flow  over  a 


bottom.  In  connection  with  the  study  of 


inverse  methods  In  section  34ot  an  explicit  example  of  a  steady 


flow  over  a  wave-shaped  bottom  was  exhibited.  However,  there 


the  surface  profile  was  given  and  the  bottom  profile  calculated. 


The  direct  problem,,  in  which  the  bottom  profile  and  other  flow 


data  are  given,  has  also  been  considered  by  several  persons. 
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Lavrent’ev  [ 1943]  announced  theorems  concerning  this  problem,  but 
did  not  include  them  in  his  later  [ 1947  ]  e x pob i i i on .  Ge  r be r 
[  1955 1  has  given  a.  comprehensive  treatment  of  the  ’’supercritical' 
case  and  has  announced  further  results  for  the  "subcritica 1"  case 
[19515],  Let  the  bottom  profile  S  be  periodic  and  symmetric 
about  vertical  lines  through  the  maxima  and  minima;  let  be 

its  intrinsic  equation  where  A  is  arc  length  measured  from  o, 
maximum  and  &  is  the  angle  between  the  tangent  and  the  x-direc- 

t  Lon.  Let  Gi  be  the  discharge  rate  for  the  fluid,  and  let  0 

u 

be  the  velocity  at  a  crest.  In  the  first  paper  he  considers 
flows  in  which  the  slope  of  the  surface  has  the  same  sign  as  that  of 

the  bottom  (we  reca  1 1  the  two  pc  - » i b ie  f lows  occurring  In  the 
linearized  theory  of  section  20:*).  Gerber  shows  that  there 
exists  at  least  one  solution  of  this  type  provided  the  following 
inequalities  are  satisfied  in  the  interval  between  a  maximum  and 
the  first  minimum  to  the  right: 


L zL  -t  wax  M 

i 1 

. jir  +  tx  £  ew  £  Oj 

<4* 


(35.3) 


where  £  and  £t  are  arbitrary  small  but  positive  quantities. 
If  certain  other  Inequalities,  further  limiting  g.  Q.  / ,  are 
satisfied,  he  is  also  able  to  prove  uniqueness  provided  0(4-}  0  , 


In  the  second  paper  he  announces  that  there  exists  at  least  ono 


solution  such  that  the  profile  has  slope  of  opposite  sign  to  that 


of  the  bottom  if 


(35.4) 
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and  Q/L-0%^  and  &  / Q,  are  small  enough;  here  L0  is  the 

arclength  from  a  maximum  to  a  minimum  of  5  and  &  is  the  ver . 

tical  distance.  Gerber's  methods  are  topological  (Schauder- 
Leray  theory)  and  do  not  yield  effective  methods  of  approximation, 
Moiseev  [1957]  has  also  considered  this  problem.  By  a  . . edi¬ 

fication  of  the  method  used  to  derive  Nekrasov's  Integral  equa¬ 
tion  (32.104),  he  derives  a  pair  of  nonlinear  integral  equations 
to  which  the  Lyapunov -Schmidt  method  is  applicable.  Let  c  be 
the  average  velocity  defined  by  (7.5)  for  an  allowable  value  of 
H,  (thus  P  increases  by  c\  over  a  wave  length),  and  let  G  be 
the  discharge  rate.  Then  Moiseev  finds  that  there  exists  a  se¬ 
quence  of  ve loci  ties  c .  >  ^  >  •  >0  associated  with  the  eigenvalues  of 
a  certain  linear  operator,  such  that,  if  C  tmc*  ,  there  exists 
a  unique  flow  provided  the  slope  of  the  bottom  is  sufficiently 
small.  Also,  if  C  »c(  or  etc.  then  the  solution  is  such 

that  the  s lopes  of  bottom  and  surface  are  of  the  same  sign;  if 
C,  c  C  c  *  the  slopes  are  of  opposite  sign. 


F low  ov e r  a_  b o 1 1; om  w i t h . a . declivity.  Let  the  flow  be  from 

left  to  right  and  suppose  the  bottom  profile  to  be  asymptotic  to 


horizontal  lines  as  x  -*  ,  the  one  on  the  right  being  lower 

than  that  on  the  left.  The  discharge  rate  Q,  and  velocity  e 


at  X-sr-oo  should  then  be  sufficient  to  determine  the  flow. 


The  existence  of  a  steady  flow  under  these  circumstances  has  been 


investigated  by  Haimovici  [1935]  and  Gerber  [1955],  The  former 
derives  a  pair  of  nonlinear  integral  equations,  similar  to 
Nekrasov's,  relating  9  and  T  of  (32.86).  An  Iterative  method 
is  used  to  prove  the  existence  of  a  solution.  Gerber  makes  use 

Ml 
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Again  of  the  Schaude r-Leray  theory,  The  theorems  established 
by  each  are  very  similar ,  but  Gerber’s  is  sharper,  Let  the 
bottom  be  g  1  von  intrinsically  by  S^(A^)  ,  measured  from  some 

i i xod  po 1 n t ,  Then  a  s o 1 u t i on  ex is ts  1  f 


% 


Q 


Cr 


/  j 

^  i 


9M 


i* 


i*',  f'X 


3  fa )  I  "4* 


I  'wr  I 


»G 

•l* 

*.N**,*H*M»NI*< 


<  [ 


A  6. 


u  4  l 


A 


*4  * 


j"*\ 

w 


(35,5) 


The  last  condition  assures  a  rapid  approach  to  the  horizontal 
asymptotes*  The  case  of  subcrttical  flow  does  not  appear  to 
have  been  treated  in  the  published  literature. 
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intensity  f  in  a  channel  of  depth  k  when  the  exact  boundary 
conditions  on  the  free  surface  are  retained.  If  c  is  the  velo¬ 
city  far  upstream  of  the  vortex,  ue  proves  existence  and  unique¬ 
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zontal  planes.  The  question  naturally  arises  as  to  whether  one 
can  establish  the  existence  of  such  waves  when  the  exact  boun¬ 
dary  conditions  at  the  interfaces  ar©  observed*  Kochin  [1927] 
extended  the  methods  of  Levi-Civita  and  Strulk  to  this  problem 
and  established  the  existence  of  (necessarily  symmetric)  inter- 
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Let  us  suppose  that  a  coord  tar.  *»  system  has  been  chosen  so 
that  we  may  treat  the  wave  motion  as  a  steady  flow  to  Mie  right 
Although  we  do  not  assume  the  motion  to  be  lrrotational,  there 
will  still  exist  a  stream  function  by  virtue  of  the 

continuity  equation.  The  vorticity  of  the  flow  will  be  given 
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by  - ,  and  wince  by  a  classical  theorem  the  vorticity  is 
constant  along  a  streamline,  the  following  equation  must  be 
satisfied  by  V  : 
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are  also  shown  to  be  symmetric  about  vertical  lines  through  crest 
or  trough.  She  also  demonstrates  that  among  this  class  of  waves 
for  finite  depth  there  is  a  unique  analogue  of  the  Gerstner  wave, 
in  the  sense  that  the  trajectories  of  individual  particles  are 
all  closed.  This  wave  has  recently  been,  investigated  by 
Kravtchenko  and  Dauber t  -1957),.  The  development  of  means  of  cal¬ 
culating  rotational  waves  has  been  the  subject  of  a  recent  In- 
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permanent  type  are  not  possible  in  a  heterogeneous  Ha  id,  but 
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